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ABSTRACT 


We present results from dynamical models of core-collapse supernovae in one spatial dimension, employing 
a newly-developed Boltzmann neutrino radiation transport algorithm, coupled to Newtonian Lagrangean hydro¬ 
dynamics and a consistent high-density nuclear equation of state. The transport method is multi-group, employs 
the Eeautrier technique, uses the tangent-ray approach to resolve angles, is implicit in time, and is second-order 
accurate in space. We focus on shock breakout and follow the dynamical evolution of the cores of 11 Mq, 15 Mq, 
and 20 Mq progenitors through collapse and the first 250 milliseconds after bounce. The shock breakout burst is 
the signal event in core-collapse evolution, is the brightest phenomenon in astrophysics, and is largely responsible 
for the initial debilitation and stagnation of the bounce shock. As such, its detection and characterization could 
test fundamental aspects of the current collapse/supernova paradigm. We examine the effects on the emergent 
neutrino spectra, light curves, and mix of species (particularly in the early post-bounce epoch) of artificial opacity 
changes, the number of energy groups, the weak magnetism/recoil corrections, nucleon-nucleon bremsstrahlung, 
neutrino-electron scattering, and the compressibility of nuclear matter. Eurthermore, we present the first high- 
resolution look at the angular distribution of the neutrino radiation field both in the semi-transparent regime and at 
large radii and explore the accuracy with which our tangent-ray method tracks the free propagation of a pulse of 
radiation in a near vacuum. Einally, we fold the emergent neutrino spectra with the efficiencies and detection pro¬ 
cesses for a selection of modern underground neutrino observatories and argue that the prompt electron-neutrino 
breakout burst from the next galactic supernova is in principle observable and usefully diagnostic of fundamental 
collapse/supernova behavior. Though we are not in this study focusing on the supernova mechanism per se, our 
simulations support the theoretical conclusion (already reached by others) that spherical (ID) supernovae do not 
explode when good physics and transport methods are employed. 

Subject headings: supernovae, neutrinos, radiative transfer, spectra 


The most dynamical phase of core-collapse supernova evo¬ 
lution is the epoch of core bounce, shock wave formation, and 
early shock propagation. Within a few milliseconds (ms) of 
birth, the shock wave propagates down the density gradient 
from where it is formed at ~ 10'"^ gcm“^ and ~10 kilome¬ 
ters to where it dies into an accretion shock at ^ 10'° gcm“^ 
and ~ 100-120 kilometers. Between birth and death, the shock 
wave, born in the neutrino-opaque region, emerges into the 
neutrino-transparent region. In the process, it “breaks out” 
through the neutrinospheres of the various neutrino species 
{Vg, Vg, and “u^”s), and, thereby, produces an ultra-luminous 
burst of neutrinos brighter than any other phenomenon in as¬ 
trophysics. This short transient is responsible for sapping the 
shock of energy and aborting the supernova explosion within 
10-20 ms of its birth and is a primary reason the supernova the- 
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1. INTRODUCTION 


orist has shifted focus to the “long-term” mechanism of shock 
revival in its various forms (Bethe & Wilson 1985; Herant et al. 
1994; Burrows, Hayes, & Eryxell 1995; Eryer et al. 1999; Eryer 
& Heger 2000; Rampp & Janka 2000; Janka 2001; Liebendor- 
fer et al. 2001; Liebendorfer, Mezzacappa, & Thielemann 
2001 ). 


Given the central importance of the breakout phenomenon 
to the dynamics of core collapse (as we think we understand 
it), and given the possibility that the next galactic supernova 
will be detected by neutrino telescopes far more sensitive than 
the underground detectors IMB and Kamioka II which caught 
SN1987A, we have made a detailed study of the breakout phe¬ 
nomenon, its neutrino signature, and its dependence on various 
parameters, progenitors, and neutrino microphysics. Neutrino 
breakout is the signature characteristic of the theory of collapse 
and supernovae, that if measured would speak volumes con¬ 
cerning the dynamics of a dense core otherwise obscured from 
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view. 

In §2, we describe our multi-group, multi-angle radia¬ 
tion/hydrodynamics code SESAME^ and the microphysics we 
use in this survey of the breakout phenomenon. We include a 
discussion of our implementation of the Lattimer-Swesty high- 
density nuclear equation of state (LSEOS) in Appendix A. In 
§3, we present our baseline core-collapse simulation using an 
11 Mq progenitor from Woosley & Weaver (1995). We follow 
the discussion of this reference model with studies of the de¬ 
pendence on spectral resolution (§4), progenitor mass (§5), the 
nuclear equation of state (§6), the weak magnetism/recoil cor¬ 
rection (§7.2), and the role of nucleon-nucleon bremsstrahlung 
(§7.3). In §7.4, using a new explicit algorithm for handling it 
(see Appendix B), we explore the effect of inelastic neutrino- 
electron scattering. In §8, event rates and signatures in the un¬ 
derground telescopes Super Kamiokande, SNO, and ICARUS 
are calculated and described. It is shown that the breakout burst 
is in principle detectable in all three of these neutrino observa¬ 
tories. Changes in the observed neutrino signal due to modifica¬ 
tions in the microphysics employed in the simulations are also 
presented. We wrap up in §9 with a summary and general con¬ 
clusions on the breakout phenomenon and its model dependen¬ 
cies. Our focus in this paper is on the breakout phenomenon, 
not the supernova itself. Hence, our simulations are restricted 
(for the most part) to the first 250 ms after bounce. 

2. NUMERICAL METHODS 

2.1. Radiation Transport 

The mass densities, temperatures, and compositions that ob¬ 
tain during core collapse conspire to produce regions in the 
deep core of a supernova where the mean-free-path for neu¬ 
trinos is just meters and the transport could be handled using 
simple diffusion theory. Simultaneously, there are regions at 
larger radius where each neutrino species decouples from the 
matter and free-streams to infinity. In between, at modest opti¬ 
cal depth, where the neutrino mean-free path is a modest frac¬ 
tion of the size of the object as a whole, the full transport prob¬ 
lem must be solved in order to obtain accurate values of the 
local neutrino flux and energy density - both critical in deter¬ 
mining the energy deposition profile and, hence, the subsequent 
evolution. 

Our Boltzmann transport solver (Burrows et al. 2000) is due 
primarily to Eastman & Pinto (1993), who developed an al¬ 
gorithm for solving the comoving transport equation using the 
Eeautrier variables and the tangent-ray method for angular dis¬ 
cretization. The transport solver is implicit in time and second- 
order accurate in space. The zeroth- and first-moment equations 
are iterated with the Boltzmann/transport equation using Ac¬ 
celerated Lambda Iteration (ALI) to speed convergence of the 
radiative transfer solution (see e.g., Mihalas & Mihalas 1984). 
Only the diagonal and the adjacent off-diagonal terms of the 
A matrix are used. This method is akin to the standard variable 
Eddington factor (VEE) approach, but no ad hoc flux limiters or 
artificial closures are necessary. The iteration scheme generally 
converges to a part in 10® in 2 to 10 steps, and automatically 
conserves energy in the transport sector. The Eeautrier scheme 
can transition to the diffusion limit seamlessly and accurately 
and the tangent-ray method automatically adapts with the La- 
grangean hydrodynamic grid as it moves. In constructing the 
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tangent rays, we cast them from every outer zone to every in¬ 
ner zone. Hence, if there are 300 radial zones, the outer zone 
has 299 angular groups in each quadrant of the unit circle. Be¬ 
cause of the spherical nature of the core-collapse problem and 
the need to accurately reproduce the angular distribution of the 
radiation field that transitions from the opaque (inner) to the 
transparent (outer) regions, such fine angular resolution is use¬ 
ful, though computationally demanding, as radiation becomes 
more and more forward-peaked. 

The coupling of the neutrino radiation to matter is done im¬ 
plicitly in operator-split fashion, using the Accelerated Lambda 
operator to obtain the derivative ‘dJ/dS’ (where S is the source 
function) after each hydrodynamic update. This derivative is re¬ 
quired for the implicit updates of temperature (T) and electron 
fraction (U) (Burrows et al. 2000). 

We solve the transport equation for neutrinos, i>e neutri¬ 
nos, and neutrinos, where stands for I'r, and 

Dr- Although in principle the neutrino microphysics is dif¬ 
ferent between i^r, and Dr we make the assumption 

standard in the field that they can be treated identically. The 
energy grouping is different for each species. Eor i/g neutri¬ 
nos, 1 < (MeV) < 320. Eor D^ neutrinos, 1 < (MeV) < 
100. Both are grouped logarithmically. Eor neutrinos, 
1 < Spij, (MeV) <100 with uniform grouping in energy. Our 
baseline models employ 40 energy groups for each species. In 
§4, we compare our results with calculations employing just 10 
and 20 energy groups. Our standard set of microphysics in¬ 
cludes all of the emission/absorption and scattering processes 
described in Burrows (2001) plus weak-magnetism and recoil 
modifications to the cross sections for neutrino-nucleon absorp¬ 
tion and scattering, as per Horowitz (1997,2002). In addition, 
we include inelastic neutrino-electron scattering as described in 
Appendix B. In our baseline models for this paper we do not in¬ 
clude inelastic neutrino-nucleon scattering, but include it only 
as a purely elastic scattering opacity source. 

2.2. Hydrodynamics 

During most epochs in supernova evolution, the radiation 
pressure term is negligible. However, during breakout when 
electron neutrino luminosities approach 10®^ erg s“' locally, this 
term can be important. Since there is as yet no natural way in 
PPM-like hydro codes to include the radiation pressure term in 
the solution to the Riemann problem, we have opted to couple 
our transport algorithm to a Newtonian and Lagrangean predic¬ 
tor/corrector hydrodynamics scheme that uses artificial viscos¬ 
ity for shock resolution. During and just after shock breakout, 
since the neutrinospheres and the regions of spectrum forma¬ 
tion are found at 30 to 100 kilometers, general-relativistic ef¬ 
fects are small to modest, GM/rc^ being just 3 to 10%. Nev¬ 
ertheless, particularly at later times, general-relativistic effects 
can introduce ~10% corrections (Bruenn, De Nisco, & Mezza- 
cappa 2001; Liebendorfer et al. 2001). 

In the core at bounce the zone size in our calculations is ~200 
meters. The sound speed at nuclear densities approaches ~ c/3. 
Hence, our Courant-limited explicit timestep is ~ 10“® seconds. 
This implies that to evolve a calculation for ~1 second would 
require ~ 10® timesteps. The solution to the transport equa¬ 
tion for all species is computationally expensive and on current 
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single-processor machines 10® timesteps is large. In order to 
circumvent this problem we sub-cycle the hydrodynamics after 
shock breakout and stall. After this most dynamical phase is 
over we take 2-4 hydro steps per transport step, thus cutting 
the total computational requirements significantly. Introducing 
sub-cycling ~20 milliseconds after bounce seems to introduce 
errors of less than 1% in all radiation and hydrodynamic quan¬ 
tities. 

For all of our simulations, we employ the tabular equation of 
state described in Appendix §A. 

3. THE BASELINE MODEL RESULTS 

We take as a reference model to establish the general con¬ 
text of the study an 11 Mq progenitor from Woosley & Weaver 
(1995). The model employs 300 mass zones out to M = 1.43 
Mq, where M is the interior mass; in this progenitor, at M 
= 1.43 Mq, r equals 2500 km. For the simulations presented 
here, which focus on collapse and the early (~ 100-200 ms) 
post-bounce evolution, we find this zoning sufficient. In §6, 
we explore the effects of changes in the nuclear compressibility 
modulus, but for our baseline model we use k = 220 MeV in the 
LSEOS (Appendix A). 

Figure 1 shows velocity v (in cm s“', upper left panel), 
log[g[p] (in g cm“^, upper right panel), entropy s (in k^ 
baryon“', lower left panel), and temperature T (in MeV, lower 
right panel), as a function of mass coordinate in Mq at five 
snapshots in time in our 11 Mq baseline model. Fig. 1 is to be 
compared with Fig. 2, which shows the electron fraction (Yg) as 
a function of mass in the same model. The thin solid line in 
Figs. 1 and 2 is the initial configuration. The thick solid line is 
at hydrodynamical bounce (approximately 200 ms after we start 
the calculation). We generally quote times relative to bounce. 
Bounce coincides with the core density reaching ~ 2- 3 x 10^"^ 
g cm“^, where the EOS stiffens dramatically. Note that E, in 
the very core reaches its minimum value before bounce. At 
core densities of ~ 2 x 10^^ g cm“^ the electron-neutrino out¬ 
ward diffusion speed becomes smaller than the the inward col¬ 
lapse speed and the neutrinos are effectively trapped. Soon af¬ 
terward, the process primarily responsible for the decrease in 
Ye, e~p nve, is balanced by its inverse and beta equilibrium 
is achieved (Bruenn 1985; Mezzacappa & Bruenn 1993abc). 
Although we have included inelastic neutrino-electron scatter¬ 
ing in our baseline model, the long dashed line in Eig. 2 shows 
the bounce Yg profile if one turns off this equilibration process. 
Because small changes in the Yg profile at bounce can affect 
the shock energetics significantly (Burrows & Lattimer 1983) 
and because inelastic neutrino-electron scattering affects the 
approach to beta equilibrium by modifying the neutrino phase 
space occupancy, inelastic neutrino-electron scattering should 
be included in a full treatment (Bruenn 1985; Mezzacappa & 
Bruenn 1993abc). We discuss this process and its effects on 
breakout more fully in §7.4. 

At bounce, a strong sound wave is formed deep in the core, 
which propagates down the density gradient set up by the in¬ 
falling outer stellar material. The sound wave steepens into a 
shock near M~0.6 Mq, just ^1 ms after bounce. The dotted 
line in Eig. 1 shows the shock fully formed and at maximum 
positive velocity (^1.9 x 10® cm s“^). Of course, the shock 
is associated with the generation of entropy and an increase in 
temperature, as seen in the lower two panels of Eig. 1. As the 


shock moves outward in mass and radius it dissociates nuclei 
into free nucleons. Although the matter is hot in this region, the 
electrons are still quite degenerate and the process e~p nvg 

builds up a sea of electron neutrinos that are trapped and ad- 
vected with the matter. Eigure 3 shows the luminosity of elec¬ 
tron neutrinos as a function of mass and time from bounce all 
the way through breakout. Although the shock does not form 
until M~0.6Mq, a phase-transition front which liberates free 
nucleons from nuclei moves in tandem with the sound wave 
generated at bounce. This causes the initial luminosity peak 
seen in Eig. 3 at M~0.5Mq, or, in this model, at a radius of 
just 11km (labeled ‘Initial rise due to phase transition’). The 
sudden drop in at M~ 0.6-0.7 Mq comes from a sharp de¬ 
crease in Xp, coincidentally nearly simultaneous with shock for¬ 
mation. Just as the shock is forming, but before the temperature 
rises dramatically, as shown in the lower right panel of Eig. 1 
(dotted line), a region of heavy nuclei exists and, due to the fact 
that the absorption cross section is not continuous across the 
phase transition, there is a decrease in at M~ 0.6-0.7 Mq. 

In Eig. 3, the specific snapshots in time between the 1 ms (dot¬ 
ted) line and the 17 ms (short dashed) line in Eigs. 1 and 2 are 
labeled. At 2.4 ms after bounce the spike reaches a local max¬ 
imum at M~ 0.85Mq, with exceeding 1.4 x 10®“* erg s“^ 
This spike roughly denotes the position of the shock. It is at 
this moment that the shock crosses the Vg neutrinosphere 
The neutrinospheric radius (Rf) for any neutrino species of en¬ 
ergy El, is set approximately by the following integral criterion: 

/.oo 2 

Tu{Ry,ey}= j Ky{£y,r)p{r)dr= (1) 

where is the total opacity. At just 3.2 ms after bounce, one 
can see clearly that the breakout pulse of electron neutrinos is 
beginning to move ahead of the shock. Having moved into a re¬ 
gion where the optical depth is below ~2/3, the neutrinos begin 
to free-stream. 

The decoupling of the radiation during breakout is perhaps 
better seen in Eig. 4, which shows the temperature profile 
(dashed line) with the Vg luminosity 3.2 ms after bounce (com¬ 
pare with Eig. 3). The four insets (thin solid lines on dotted 
axes) show polar plots of the specific intensity, for = 

12 MeV. The large dots on the temperature profile show the po¬ 
sition, in units of enclosed mass, at which the specific intensity 
has been sampled. At point ‘A’ (M ~0.674 Mq, r ~ 21.8 km), 
lufO) is virtually a circle, indicating that the radiation field is 
isotropic, that the net flux is nearly zero, and that the neutrinos 
at that position and energy are trapped. In contrast, at point ‘D’ 
(M ~1.03 Mq, r ~ 124.0km), is very forward-peaked 
and well-resolved by our tangent-ray algorithm with 160 an¬ 
gles in each quadrant of the unit circle. Insets ‘B’ and ‘C’ show 
lufO) at intermediate positions, on either side of the shock. The 
small increases in in the forward quadrants of ‘C’ are 

consistent with limb-brightening effects from the shock, never 
before resolved. 

Eigure 3 shows a marked decrease in the breakout Vg lumi¬ 
nosity pulse just as it begins to propagate ahead of the shock, 
between the lines labeled ‘t=2.4ms’ and ‘t=6.8ms’. The lo¬ 
cal maximum of at 2.4 ms occurs at a radius of approxi¬ 
mately Tpeak — 41 km. We expect the peak luminosity to 

be attenuated approximately by an amount ~ exp(-Ta) in prop¬ 
agating to infinity (L“), where is the total absorptive opti¬ 
cal depth exterior to rpeak- >Erom the results presented here, 
we find Ta(rpeak) ~ 1.9 at ~ 12 MeV. This simple estimate 
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shows that we should expect to decay by ~ 85%. With 
~ 1.45 X 10^“* erg s“' we expect ~ 0.22 x 10^'^ erg s“' 
and find from the simulation that L“ ^ 0.25 x 10^“* erg s“*, in 
good agreement with our estimate. This asymptotic luminosity 
is shown clearly as the bump in the last line in Fig. 3, labeled 
‘Breakout Pulse’, which has moved out to 1.35 Mq. Behind the 
pulse, at M~ 1.08 Mq , is a sharper peak in the luminosity that 
denotes the shock position. During this breakout phase of shock 
evolution, much of the shock energy is sapped by the neutrino 
losses that attend electron capture on newly-liberated free pro¬ 
tons. This process also creates a characteristic deleptonization 
trough manifest by the marked decrease in Ye in Fig. 2 near 
M~0.8Mq. The shock stalls quickly and in just a few mil¬ 
liseconds after bounce all velocities are < 0. Note the negative 
entropy gradient between M~ 0.75 Mq and M~ 0.95 Mq. De¬ 
pending upon the compositional gradients in this region, in a 
multi-dimensional simulation this region might be convectively 
unstable (Burrows and Fryxell 1992; Keil, Janka, and Muller 
1996). 

Although the shock continues to move outward in mass as 
matter flows through it from the free-falling outer core, all 
matter velocities are negative by the time the shock reaches 
r ^ 80-90km. The short dashed lines in Figs. 1 and 2 show 
the configuration of the model 17 ms after bounce and ~10ms 
after the shock has stalled. The dot-dashed line shows these ba¬ 
sic hydrodynamical quantities 200 ms post-bounce. The shock 
has consumed another 0.2 Mq of infalling material and the en¬ 
tropy behind the shock has increased dramatically. The tem¬ 
perature over the whole post-shock profile has increased as a 
result of the compressional work done on these zones by the in¬ 
falling stellar material. In particular, at the peak in the temper¬ 
ature profile at M^0.7Mq, T increased by ^4 MeV between 
f = 17 ms and t = 200 ms. In these hundreds of milliseconds 
after bounce, the shock moves outward in radius to ~ 175 km. 
Although it is in this epoch that neutrino heating behind the 
shock is thought to revive it, no explosion is seen in our sim¬ 
ulations at this time. Figure 5 shows the evolution of selected 
mass zones in our 11 Mq model as a function of time. Small 
oscillations are visible just after bounce - the ‘ringing’ of the 
core - and the shock position is made clear by the change in 
the infalling matter velocity between 100 km and 200 km dur¬ 
ing the 280 ms of post-bounce evolution shown here. Although 
the neutrino-driven mechanism is not the focus of this work, 
from this figure and from many similar calculations we see no 
evidence of explosion. 

In Fig. 3, we showed as a function of mass at various 
snapshots in time. The last snapshot shows the Ve breakout 
pulse at 1.35Mq and propagating at the speed of light to the 
edge of the grid. Complementary to Fig. 3, the time evolution of 
at infinity is shown in Fig. 6 (thick solid line). The luminos¬ 
ity of anti-electron neutrinos (Lp^, thin solid line) and the com¬ 
bined luminosity of and Vr (collectively, dotted 

line) are also shown for the first 250 ms of post-bounce evolu¬ 
tion. Note that just before the major breakout pulse of electron 
neutrinos, there is another small peak that reaches ~ 0.72 x 10^3 
erg s“'. This peak results from the deleptonization of the core as 
the Ye drops in Fig. 2 from that of the initial model to ij. ~ 0.27. 
There is a small dip in just after this initial rise. This dip in 
the asymptotic luminosity signals the end of the infalFcollapse 
phase of the supernova. After the initial increase in Lp^ due to 
deleptonization, Lp^ decreases because of the increased opacity 


in the dense collapsing core. The high opacity decreases the lo¬ 
cal luminosity at small radii and isolates the precursor peak in 
luminosity at larger radii caused by early deleptonization. If it 
were not for shock formation, Lp^ at infinity would continue to 
decline. The main breakout pulse just after the small downturn 
in Lp^ results from the dissociation of nuclei and subsequent 
electron capture on free protons, as described above and shown 
in Fig. 3. 


In Fig. 7 we show the corresponding average neutrino energy 
for each species over the same time interval as Fig. 6. There 
are many ways one might define the average neutrino energy. 
The literature does not consistently employ a unique prescrip¬ 
tion and considerable confusion exists. The average energy we 
present in our figures is the rms average energy, defined by: 

■/ <i£p47p(£p)1 


= 


( 2 ) 


J" dSyJyi^Sy) j 

where JuiSu) is the zeroth moment of the specific intensity. We 
use this definition in order to make comparisons with other re¬ 
cent work (Liebendorfer et al. 2001; Rampp et al. 2002). After 
approximately 50 ms, the average energies show a nearly linear 
increase with time. Mezzacappa et al. (2001) and Liebendor¬ 
fer et al. (2001) find similar behavior. In addition, although we 
employ a different progenitor, the values for the average neu¬ 
trino energies they obtain compare well with those presented 
in Fig. 7. For example, 250 ms after bounce in their 13Mq 
progenitor, Mezzacappa et al. (2001) obtain (£p ) ~ 21 MeV, 
(£pj ~ 16 MeV, and (£pj ~ 14 MeV. We find that (£p^) ~ 20 
MeV, (£p,) ~ 15.5 MeV, and (£pj - 13 MeV. 


Figure 8 shows the luminosity spectrum of electron neutri¬ 
nos at infinity in units of 10^^ erg s“^ MeV“'. The thick solid 
line denotes the maximum peak breakout pulse luminosity and 
corresponds to the 2.55 x 10^^ erg s“' peak in in Fig. 6. The 
thin solid lines in this plot show the rise up to the peak and 
the simultaneous shift in {e^^) depicted in Fig. 7. The dashed 
lines show the evolution of the luminosity spectrum after the 
breakout pulse. The pre-breakout lines show the spectrum at 
11.6, 5.1, and 1.3 ms before the peak in the breakout pulse. The 
dashed, post-breakout, lines show the spectrum at 4.2, 9.6, and 
40.5 ms after the Ve breakout pulse. Although in this figure the 
peak energy is dropping as the luminosity decays, the high en¬ 
ergy tail becomes broader and with our definition of the aver¬ 
age energy (eq. 2) {e^^) increases. Figure 9 shows the emergent 
spectrum much later (at 210 ms after bounce) and for all three 
neutrino species. The thick solid line shows the Ve spectrum, 
the thin solid line shows the Ve spectrum, and the dotted line 
shows the combined spectrum. The post-bounce energy hi¬ 
erarchy one expects from simple considerations of the opacity 
sources for the various species, {e^^) < (sp^) < {Sp^), is evident 
in both Figs. 7 and 9. 


Although this exact progenitor model has not been used in 
any of the recent supernova modeling done by groups employ¬ 
ing sophisticated neutrino transport (Rampp 2000; Rampp & 
Janka 2000, 2002; Liebendorfer et al. 2001; Liebendorfer, 
Mezzacappa, & Thielemann 2001; Mezzacappa et al. 2001), 
we find that the basic structures and systematics are similar 
between groups. In particular, our core entropy (Score) drops 
slightly at the beginning of the calculation and then climbs 
during bounce to Score — 1-32 kg baryon“'. We observe sim¬ 
ilar behavior for our 15 Mq progenitor with Score — 1-33 kg 
baryon“'. Without inelastic neutrino-electron scattering we ob¬ 
tain Score — 1-16 ks baryon“'. The absolute magnitude of these 
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results and the effect of neutrino-electron scattering on Score are 
duplicated in the work of Rampp (2000), whose Figs. 5.3a and 
5.6b shows that with and without neutrino-electron scattering in 
their 15M0 model they obtain Score — 1-3 and 1.14 baryon“', 
respectively. A recent collapse calculation by Liebendorfer et 
al. (2002b) of a 13 Mq progenitor also shows Score ~ 1-3 kg 
baryon"'. Core (F”"^®) in our IIMq model drops to 0.266 
and 0.293 at bounce in our models with and without neutrino- 
electron scattering, respectively. Both numbers from Rampp 
(2000) are slightly higher; his Fig. 5.3b gives the correspond¬ 
ing numbers as approximately 0.28 and 0.30-0.31, respec¬ 
tively. Liebendorfer et al. (2002b) have ~ 0.30, including 
neutrino-electron scattering. Peak positive velocities in all our 
calculations reach ~ 2 x 10® cm s“' at M~ 0.7Mq - in good 
agreement with previous work. 

One important difference bears mention. Our peak breakout 
electron-neutrino luminosity is slightly lower than that obtained 
in other work. Both Liebendorfer et al. (2001) and Rampp & 
Janka (2000) obtain of ~ 3.5 x 10^^ erg s“^ with compa¬ 
rable progenitors. Indeed, Liebendorfer et al. (2002a) present 
the breakout pulse for several different progenitors and find 
virtually identical for all models. We obtain ~ 

2.55 X 10^^ erg s ^ in our IIMq and virtually the same num¬ 
ber with the 15 Mq and 20 Mq models we present in §5. Our 
results consistently give lower by about 30%. This lu¬ 
minosity pulse is a signature of the most dynamical phase in 
supernova modeling and it is important to achieve adequate an¬ 
gular, energy, and spatial zoning to resolve it. For comparison, 
Mezzacappa et al. (2001) and Liebendorfer et al. (2001) typ¬ 
ically employ eight angular bins using the S„ method at each 
radial point and twelve energy groups. In particular, poor angu¬ 
lar resolution at large radii may compromise the propagation of 
the breakout pulse. Liebendorfer et al. (2002b) state that they 
experience significant “numerical diffusion” of their breakout 
pulse as it propagates from the point of decoupling to the edge 
of their computational grid. The peak asymptotic luminosity 
they quote, ~3.5xl0^^ erg s“\ is sampled at 500km in 
their calculations. This peak drops ~30% to ~2.25xl0^^ erg 
s“^ in propagating to ~4200km (see their Fig. 14 and Fig. 20). 
Plotting the radial coordinate of the peak in the luminosity pulse 
as a function of r-ct, Liebendorfer et al. (2002b) find that 
decreases, shifts by ~ 500- 1000 km, and broadens due to nu¬ 
merical diffusion. Constructing a similar plot with our collapse 
calculations, over a similar range of radii, we find a shift in 
the pulse peak of just <10 km, or a fractional deviation from c 
of ~0.3%. We also find a decrease in between 500 km 
and 4200km of 8.5%. In short, we obtain a lower asymptotic 
Vg breakout pulse than Liebendorfer et al. (2002b), that cannot 
be explained by the considerable numerical diffusion in their 
first-order upwind differencing scheme. Using the tangent-ray 
method, ourLjJ^*^ does not experience the same degradation and 
maintains an on-grid speed remarkably close to c after fully de¬ 
coupling from the matter. 

Although they describe it as diffusion, the temporal disper¬ 
sion and decay Liebendorfer et al. (2002b) experience in the 
electron-neutrino breakout pulse might be better interpreted as 
a consequence of the progressive loss of angular resolution as 
the luminosity pulse propagates outward and, thereby, becomes 
more and more forward-peaked. Though a fixed angular grid 
with eight points can not resolve a sharply-pointed pulse, it 


should be adequate for calculating the flux factors of relevance 
to neutrino energy deposition in the gain region and for gaug¬ 
ing the viability of the neutrino mechanism of supernova explo¬ 
sions. 

Two hundred milliseconds after bounce, the shock has moved 
out to about 200km in radius. A gain region of net heating 
has formed between the gain radius Rg ~ 100 km (where heat¬ 
ing balances cooling) and the shock. For = 10.7 MeV, the Vg 
neutrinosphere {Ru,, as defined in eq. 1) occurs at ~ 60 km. Be¬ 
tween R^^ and Rg is a region of net cooling. The solution to the 
Boltzmann equation yields the specific intensity of the neutrino 
radiation field, /j^, as a function of both energy and angle, at ev¬ 
ery radial coordinate and timestep. When neutrinos of a given 
energy are trapped and diffusive, their radiation field is isotropic 
in angle. Beyond the neutrinosphere, where they approach 
the free-streaming limit, the neutrino radiation field becomes 
continuously more forward-peaked. In an effort to understand 
the angular distribution of Fig. 10 shows a polar plot of 
/i^^(0) for many Vg energies at R^S^v^ = 10.7MeV) = 60km. All 
curves in this figure have been normalized to the same maxi¬ 
mum specific intensity. The outer near-circle (thick solid line) 
in this figure shows for = 320 MeV. Fig. 10 shows 
that although 60km is R^^ for neutrinos with = 10.7 MeV, 
320MeV neutrinos are completely trapped at this radius; /i/,(6*) 
is isotropic and the net flux is nearly zero for = 320 MeV. 
This is a consequence of the fact that the neutrino opacity is 
proportional to e^. Also shown is 1^(9) for other neutrino en¬ 
ergies down to 1 MeV (also a thick solid line). Only every 
other neutrino energy calculated is shown. Going from high 
to low energies, one sees that the radiation field becomes less 
isotropic, the lowest energies being the most forward-peaked. 
The = 10.7 MeV point is highlighted (middle thick solid 
line) to show the angular distribution of the radiation field at 
the neutrinosphere for that energy. A glance at Fig. 10 shows 
how forward-peaked the radiation field is at the point of de¬ 
coupling (the neutrinosphere), for a given energy (in this case, 
10.7 MeV). Note that in this calculation, R^^ occurs at our 270th 
radial zone, so that in each quadrant of the polar plot there are 
269 angular bins. 

Also of interest is the dependence of the angular distribu¬ 
tion of the radiation field on radius for a given energy. Figure 
11 shows how = 10.7 MeV) evolves from the neutri¬ 

nosphere radius (~60 km, outer thick solid line) to the shock ra¬ 
dius (Rs ~ 200km, inner thick solid line). Also shown is 1^X9) 
at Rg, where heating balances cooling, and /i/,(6*) at larger radii, 
where heating dominates. This figure shows that the angular 
distribution of the radiation field is quite forward-peaked in the 
heating region. As noted in Liebendorfer et al. (2002b), insuf¬ 
ficient angular resolution in this region can lead to artificial en¬ 
hancements in the energy deposition profile and, hence, could 
affect the subsequent dynamics. In Fig. 12, we show the neu¬ 
trinosphere radius (as defined in eq. 1) as a function of energy 
for all neutrino species at four different snapshots in time from 
just after breakout to 200 ms after bounce. R^^,Rii^, and R,^^ are 
the thick, thin, and dotted lines, respectively. The break in the 
neutrinosphere radius profiles at approximately 50 km, 70 km, 
and 100 km in the first three panels roughly denotes the posi¬ 
tion of the shock. The bottom-right panel corresponds roughly 
in time with the polar plots of the radiation field in Figs. 10 and 
11. Note that the energy hierarchy, < (Ej/^), is 

reflected here; R^^ < < R^^. 
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4. SPECTRAL RESOLUTION 

Our baseline model employs 40 energy groups for each neu¬ 
trino species, with the actual grouping described in §3. In an 
effort to understand the effects of degrading the spectral resolu¬ 
tion, we have run the same 11 Mq model, with the same spatial 
zoning and the same physics as in the baseline model, but with 
20 and 10 energy groups. All three calculations were carried 
out to 250 milliseconds after bounce. 

We observe only small quantitative differences and virtually 
no qualitative differences in the overall evolution of all quan¬ 
tities in comparing the run with 40 energy groups to that with 
only 20. Over the whole post-bounce evolution, is ~0.05 
MeV higher in the 20-group calculation, a difference of just 
0.5% at 250 ms. However, is higher by 2% in the model 
with just 20 energy groups. Differences in the average energy 
and luminosity of Vg and neutrinos amount to less than ~ 1% 
overall, with each of these quantites lower in the 20-group cal¬ 
culation than in the model run with 40 energy groups. The basic 
hydrodynamical quantities were similarly unaffected. The core 
entropy was reproduced in the the 20-group model to within 
0.5%, with the higher resolution model having slightly lower 
entropy. The core Ij, at bounce was just 0.5% lower in the 
model with 20 groups. The only qualitative difference observed 
between the two models was a slight oscillation in the velocity 
versus mass profile inside 0.5 M© just before bounce. 

We find significant qualitative and quantitative differences 
between our 10-group model and our models with higher spec¬ 
tral resolution. The small oscillations observed in the matter 
velocity in the 20-group model are much more pronounced on 
infall with just 10 groups. There are similar and seemingly cor¬ 
related oscillations in Ij,. Oscillations very similar to these have 
been observed in Rampp (2000) and Mezzacappa & Bruenn 
(1993c). The oscillations result from insufficient spectral res¬ 
olution of the Fermi surface of electron neutrinos. In addition 
to the oscillations, many components of the hydrodynamic and 
neutrino spectral evolution are affected by the use of only 10 
energy groups in our calculation. At bounce, the core entropy, 
which had come within 5% of that obtained in the 20- and 40- 
group calculations (Score ~ 1 -35 ks baryon”'), dropped to 0.8 kg 
baryon“*. This, of course, also manifests itself in the temper¬ 
ature profile. The core temperature after bounce reaches only 
7.5 MeV, as compared with the 11 MeV obtained in the cal¬ 
culations with 20 and 40 energy groups. Despite the fact that 
the core Ij, was slightly higher in the 10-group calculation, the 
peak positive shock velocity just after bounce was 15% lower. 
The Ve breakout pulse was 10% lower in the 10-group model 
than in the 40- or 20-group models. was ~ 0.03 x 10^^ 
erg s“' higher in the low-resolution model throughout the post¬ 
bounce evolution, this difference amounting to ~ 8% at 250 
ms. In addition, we find that Lp^(0 in the 10-group model peaks 
faster and higher than the 40-group model, but drops off more 
quickly in time. Although at 250 ms is just 9% lower in 
the 10-group model than the 40-group model, the slope of the 
luminosity is steeper and that difference grows with time. The 
average energy for all species is higher in the 10-group calcu¬ 
lation than in the 40-group calculation, but not by more than a 
few percent over the whole post-bounce evolution. 

In sum, we conclude that 10 neutrino energy groups are in¬ 
sufficient for resolving the neutrino radiation field. We see the 
largest differences in the emergent spectrum for the electron 


neutrinos. This is likely due to the fact that the grouping 
must extend to ~ 300 MeV in order to resolve the Fermi energy 
of the electron neutrinos within the core. Because Ve transport, 
through the charged-current interaction, affects the dynamics so 
significantly, failure to resolve the Ve radiation field affects most 
observables. Our model employing 20 energy groups comes 
quite close to reproducing our results with 40 energy groups. 
Differences are no larger than 2% and, for most quantities, the 
fractional differences are closer to 0.5-1%. Considering the 
large decrease in computation time (a factor of 2), 20 groups 
are probably desirable for all but the most detailed models. In 
comparison, with 21 energy groups Rampp (2000) finds that 
the Ve Fermi surface is under-resolved, as the large oscillations 
he sees in the trapped lepton fraction demonstrate. Typical re¬ 
cent calculations by Rampp & Janka (2002) employ 20-30 en¬ 
ergy groups spaced geometrically. Mezzacappa et al. (2001) 
and Liebendorfer et al. (2001) employ 12 energy groups. 

5. COMPARISON BETWEEN PROGENITORS: 

15 M0 AND 20 M0 MODELS 

A complete theory of core-collapse supernovae must hope to 
understand the detailed dynamics and neutrino signatures of all 
possible supernova progenitors. If successful neutrino-driven 
supernova models were obtained, one might hope to find sys¬ 
tematic trends between, for example, the progenitor mass and 
the peak neutrino luminosity, final supernova energy, or mass 
of the nascent protoneutron star. We begin such an investiga¬ 
tion by first considering the dynamics and emergent neutrino 
characteristics of two more massive progenitors. 

The top two panels of Fig. 13 show the (10^“* erg s“', up¬ 
per left-hand panel) and and (upper right-hand panel), 
at infinity as a function of time for three different progenitor 
masses: IIM 0 (the baseline model, thick solid line), 15 M© 
(thin solid line), and 20 M 0 (dotted line). Note that the Ve 
breakout pulse for each progenitor is remarkably similar, a 
result also recently obtained in the work of Liebendorfer et 
al. (2002a). Interestingly, as discussed in §3, our Vg breakout 
pulse reaches ~ 2.55 x 10^3 erg s ', whereas Liebendorfer et 
al. (2002a) obtain a characteristic peak Lp^ of ^ 3.5 x 10^3 erg 
s“'. Note that our 20 M© model, which also employs 300 spa¬ 
tial zones, 40 neutrino energy groups, and the tangent-ray al¬ 
gorithm for angular binning, reaches a somewhat lower peak 
luminosity at ~ 2.225 x 10^3 erg s“'. 

Lp^ and Lp^ are virtually identical in both the 11 M© and 
15 M 0 models over the first 100 ms of post-bounce evolution 
and only later develop small differences. In marked contrast, 
the 20 M 0 model exhibits higher neutrino luminosities for all 
species after electron-neutrino breakout. For Lp^ and Lp^, the 
difference between the 11 M 0 model and the 20 M© model is 
a factor of two at 200 milliseconds after bounce. For Lp^, the 
difference is a factor of 1.5. In addition, Lp^ and Lp^ peak ap¬ 
proximately 100 ms later in the 20 M 0 model than in the 11M© 
progenitor. 

Also shown in Fig. 13 (lower two panels) are the correspond¬ 
ing average neutrino energies (defined by eq. 2) of the emergent 
spectra. The lower left-hand panel shows the time evolution of 
(epj and the lower right-hand panel shows the same for (epj 
and (eTp^). One hundred milliseconds after bounce the differ¬ 
ences between the spectra from each model become significant. 
Not only are the spectra for the 20 M 0 model harder, but the 
average energies are increasing faster than those for the 11M© 
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and 15 Mq models. Some of these systematics can be under¬ 
stood by inspecting Fig. 14, which shows the mass flux (M) as 
a function of radius in units of Mq s“* for the 11 Mq (thick 
solid line), 15 Mq (thin solid line), and 20 M© (dotted line) 
models. Also shown are the positions of the neutrinosphere 
for the average emerging neutrino energy. Open triangles cor¬ 
respond to whereas filled squares and open circles denote 

and R^^, respectively. In the steady-state, the total luminos¬ 
ity of a given neutrino species should be given approximately 
by the accretion luminosity, written in terms of M and R^: 

~ GMyM/R^, where is the mass enclosed by Ry. This 
reflects the role of the neutrino luminosity in carrying away 
the energy of infall. To zeroth order, the higher mass flux in 
the 20Mq model in Fig. 14 explains the higher luminosity at 
that epoch in the top two panels of Fig. 13. This conclusion is 
modified slightly by the actual radial positions of the various 
neutrinospheres for each progenitor in this epoch and the mass 
enclosed within R,y. Given the simple expression for the 
ratio of the total neutrino luminosity of the 20 Mq model to that 
of the 11 Mq model should be approximately 

^2OM0 ^ /^20Mo\ /^2OM0\ /^HMqX 

We find that the actual ratio, as obtained by solving the full 
problem, is ~ 1.9. These significantly higher neutrino lumi¬ 
nosities in the post-bounce phase and slightly harder average 
energies affect the number of neutrino events we expect from 
the next galactic supernova. In §8 we present detection results 
for each progenitor. 

Although we obtain higher average neutrino luminosities and 
energies for the 20 Mq progenitor, there is no sign of explo¬ 
sion or developing explosion in this model. In fact, approx¬ 
imately 200 ms after bounce the shock in this progenitor has 
reached r ^ 190 km and begins to recede; at approximately 
250 ms, it has dropped to r ~ 150 km. In contrast, both the 
IIMq and 15 Mq progenitors maintain shock radii of 180- 
190 km throughout this epoch. Although none of these simu¬ 
lations yields an explosion, the lower mass models seem more 
promising. 

6. DEPENDENCE ON THE NUCLEAR EQUATION OF STATE 

One might expect that large modifications to the nuclear 
equation of state could have a significant impact on the collapse 
and bounce dynamics, as well as on the emergent neutrino spec¬ 
tra. We have not yet studied nuclear equations of state based 
on finite-temperature mean-field theory, nor have we incorpo¬ 
rated exotic particle species and phase transitions. Such equa¬ 
tions of state are of considerable interest, but are not always 
easily incorporated into existing dynamical codes. In addition, 
the neutrino microphysics would have to be made consistent 
with the presence of exotic constituents (e.g., neutrino-meson 
and neutrino-quark interactions) (Reddy, Prakash, & Lattimer 
1998). 

Performing an albeit limited exploration of the parameter 
space available to high-density equations of state, we have con¬ 
structed three tabular versions of the LSEOS (see Appendix § A) 
with different nuclear compressibilities: n = 180, 220, and 375 
MeV. For all of our baseline models we took k = 220 MeV. 
Note that, since Lattimer and Swesty (1991) attempted to con¬ 
struct a thermodynamically-consistentEOS with self-consistent 


Maxwell constructions at the phase transition to nuclear matter, 
LSEOS models with different compressibilities also have dif¬ 
ferent /ts, compositions, and energies, not just pressures. 

In order to illustrate the role of the “nuclear compressibil¬ 
ity," we have run collapse and bounce simulations with all three 
KS, evolving each to 200 ms after bounce. Eigure 16 shows 
the temperature profile in an 1 IMq progenitor for calculations 
with K = 180 (dashed lines) and k = 375 (solid lines) ~2ms 
and ~ 100 ms after bounce. The model with k = 180 reaches a 
higher core temperature, but has a steeper negative temperature 
gradient, reaching a minimum at M~ 0.4 Mq. The peak tem¬ 
perature at M~ 0.65 Mq is also higher for the lower compress¬ 
ibility model. Overall, however, the two models are remarkably 
similar. Peak velocities in both calculations reach ~2 x 10^ cm 
s“^ at M~ 0.7 Mq. In addition, the entropy and Yg profiles are 
qualitatively identical for all three ks. Although the i/g breakout 
luminosity pulse for the model with k = 375 is delayed by 1 ms 
compared with the corresponding pulse for both «: = 180 and 
K = 220, the shape and magnitude of the pulses are very similar. 

Eigure 17 shows that for both EOS models the luminosities 
at the edge of the grid for all neutrino species (apart from the 1 
ms offset for k = 375) evolve similarly during the first ~ 30 ms 
after bounce. The average neutrino energies for both EOS mod¬ 
els are also similar during this earlier phase. After 50 ms, the 
differences between the models with k = 180 and k = 375 in 
both the neutrino luminosities and spectra become discernible. 
The post-bounce neutrino luminosity for all species is lower 
in the k = 375 model. At 200 ms after bounce the fractional 
difference is approximately 4%, 5%, and 9% for and 

respectively. The model with k = 375 also produces char¬ 
acteristically softer spectra with the average energy for each 
species being ^3% lower. Although and (epj are lower 
by a few percent in the model with k = 375, ) is not evolv¬ 

ing to higher energies as quickly as in the model with k = 180. 
Eor all models, (sp^) evolves approximately linearly with time, 
but the model with k = 180 (the softer EOS) evolves with a 
larger slope. In §8.4 we quantify the importance of the nuclear 
compressibility modulus in shaping the detected neutrino signal 
from supernovae. 

These differences in spectra and some of the basic hydro- 
dynamical quantities are not large enough to affect the post¬ 
bounce dynamics significantly on 100 ms timescales. In fact, 
200 ms after bounce, there is very little difference between 
the shock position, gain radius, or neutrinosphere positions in 
these two models. Perhaps larger differences between mod¬ 
els with different compressibilities develop on one- and ten- 
second timescales in the post-explosion protoneutron star cool¬ 
ing epoch. This might be important both for detectability of 
the neutrino signature at late times and the evolution of the 
neutrino-driven protoneutron star wind (Thompson, Burrows, 
and Meyer 2001). 


7. MODIFICATIONS TO THE STANDARD CROSS SECTIONS 
AND SOURCES 

The work of Burrows & Sawyer (1998,1999), Raffelt & 
Seckel (1998), Janka et al. (1996), and Reddy, Prakash, & Lat¬ 
timer (1998) all point to a reduction at high density of the dom¬ 
inant neutral- and charged-current opacity sources for neutri¬ 
nos due to nucleon-nucleon correlations caused by the strong 



interaction and Fermi statistics. Such modifications to the stan¬ 
dard neutrino opacities depend on the model for strong interac¬ 
tions employed. Ideally, they should be included consistently in 
constructing the nuclear equation of state (Reddy, Prakash, & 
Lattimer 1998). Although they did not take into account these 
changes to the nuclear EOS, Rampp et al. (2002) have recently 
investigated the role of nuclear correlations in their dynamical 
models. They found modest enhancements in and and 
considerable enhancements in (~20%). Although inclusion 
of these processes did not lead to explosions, they did affect the 
dynamics and it seems they should not be absent from a full 
treatment of the problem. 

In a future work we will incorporate the dynamical struc¬ 
ture function formalism of Burrows & Sawyer (1998,1999) and 
Reddy, Prakash, & Lattimer (1998). For now, we provide a set 
of simple tests, varying some of our cross sections and reaction 
rates in order to test the sensitivity of the observables and the 
dynamics. 

7.1. Artificial Opacity Reduction 

In Liebendorfer (2000), erroneous explosions were obtained 
in one-dimensional models of supernovae because of an artifi¬ 
cially decreased neutral-current neutrino-nucleon opacity. As 
a simple test, we have decreased the neutrino-neutron ((t„) and 
neutrino-proton {ap) scattering cross sections by a factor of 10 
everywhere, for all neutrino energies. Despite such a change, 
because the dynamics during collapse and bounce are domi¬ 
nated by neutral-current scattering on nuclei and the charged- 
current process VgU ^ e~p, the breakout phenomenon is only 
modestly affected; the peak Vg breakout pulse is just 10% higher 
than in our baseline model and reaches 2.85 x 10^^ erg s“'. Two 
hundred and fifty milliseconds after breakout, ~ 0.25 x 10^^ 
erg s“' in our baseline model. In our model with (t„^^/ 10, we 
obtain ~ 0.35 x 10^3 erg s Figure 15 shows Lp^ and 
Lp^ at infinity as a function of time for the baseline model 
and the model presented here. Both Lp^ and Lp^ are increased 
by ^40-50% in the calculation with lower opacity. These 
changes, however, are modest compared with those for Lp^ and 
(£p^). Two hundred and fifty milliseconds after bounce, Lp^ 
equals 2.4 x 10^^ erg s“', a factor of four larger than in the 
baseline model. For Vp neutrinos, there is no charged-current 
absorption process on free nucleons to contribute to the total 
opacity, as there is for the Vg and Vg neutrinos. For the the 
neutral-current scattering dominates the opacity and so chang¬ 
ing that opacity source so dramatically has significant conse¬ 
quences forLp^ and (Cp^). 

Importantly, because of the increase in the Lp^ and Lp^ we ob¬ 
tain a larger gain region in this calculation, higher entropy be¬ 
hind the shock, and a larger shock radius 250 ms after bounce. 
In the baseline model, at 250 ms post-bounce the shock is sit¬ 
ting at 180 km, whereas in the model with reduced opacities, 
Rs ~ 210 km. Still, although the dynamics are affected by this 
drastic decrease in neutrino opacity, we do not obtain an explo¬ 
sion in this first 250 ms of post-bounce evolution. 

7.2. The Weak Magnetism/Recoil Correction 

Our baseline 11 Mq model includes the weak mag¬ 
netism/recoil correction to the charged-current opacities for 
Vg and Vg and to the neutral-current scattering opacities off 


of nucleons for all neutrino species (Vogel 1984; Horowitz 
1997,2002). The weak magnetism/recoil correction lowers 
the cross section for all processes; the decrease is largest for 
anti-neutrinos {vg and Vp). We ran our baseline model (§3) 
without these corrections and found only small differences in 
the first 250 milliseconds of post-bounce evolution. Lp^, Lp^, 
(£pj,and (£p^) are all lower at the 0.5% level without the weak- 
magnetism correction. (£pj and Lp^ are less by 4% in this rep¬ 
resentative 11 Mq calculation. This slight softening and dim¬ 
ming of the emergent spectra can be understood simply as the 
result of increasing the total opacity for all species by removing 
the small weak magnetism/recoil correction. Hence, though the 
corrections themselves for given temperatures, densities, and 
neutrino energies can be greater than 10-15%, due to feedbacks 
on the dynamics and thermodynamics, the self-consistent inclu¬ 
sion of weak-magnetism/recoil effects in a core-collapse calcu¬ 
lation has only very limited consequences. 

7.3. The Role of Nucleon-Nucleon Bremsstrahlung 

The effects of bremsstrahlung on the Vg and Vg spectra are 
negligible at virtually all times because for them the charged- 
current processes Vgn ^ pe~ and Vgp ^ ne^ dominate. Thomp¬ 
son, Burrows, and Horvath (2000) and Raffelt (2001) showed 
that nucleon-nucleon bremsstrahlung might have important 
consequences for the emergent Vp spectra and that it domi¬ 
nates e'^e~ VpVp as a production process. Figure 18 shows 
the effect of bremsstrahlung on the average energy of VpS. The 
plot itself shows the average energy (in MeV) as a function of 
radius approximately 230 ms after bounce. The inset shows 
Rp^(£p^), the neutrinosphere radius, defined by eq. (1), as a 
function of energy. The solid line is the baseline model, with 
neutrino-electron scattering redistribution and nucleon-nucleon 
bremsstrahlung. The short dashed line is the same model 
run through the full evolution, but without nucleon-nucleon 
bremsstrahlung. The long dashed line shows the same model 
without neutrino-electron scattering energy redistribution, but 
with nucleon-nucleon bremsstrahlung. One can see immedi¬ 
ately from the inset that at low energies, bremsstrahlung in¬ 
creases the opacity significantly. It is also clear from the main 
figure that inelastic neutrino-electron scattering profoundly ef¬ 
fects the emergent Vp spectrum. 

Figure 19 shows Lp^ and (£p^) at infinity as a function 
of time for the baseline model (thick solid line), which in¬ 
cludes nucleon-nucleon bremsstrahlung, and the same model 
computed without bremsstrahlung (thin solid line). With 
bremsstrahlung, the average energies of the Vp neutrinos are 
lower and the luminosity is significantly higher. The results 
of Thompson, Burrows, and Horvath (2000) anticipated these 
results based on comparisons between the bremsstrahlung and 
e'^e“-annihilation neutrino production rates. Near the Vp neu¬ 
trinosphere, the bremsstrahlung emission spectrum is peaked 
at low neutrino energies (< 10 MeV) compared with an¬ 
nihilation. For this reason, Thompson, Burrows, and Horvath 
(2000), Burrows et al. (2000), and Raffelt (2001) concluded 
that the emergent Vp spectra have both higher luminosity and 
lower average energy. As evidenced by Fig. 19, this prediction 
is borne out in our dynamical models. 

7.4. The Ejfects of Inelastic Neutrino-Electron Scattering 
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Figure 20 shows the luminosity at infinity as a function of 
time for all species in the 11 Mq models with (thin solid line) 
and without (thick solid line) inelastic neutrino-electron scat¬ 
tering, using the prescription described in detail in Appendix 
B.l. Figure 21 shows the rms neutrino energy for each species 
for the model without inelastic scattering and should be 
compared with Fig. 7. The main breakout pulse is en¬ 
hanced slightly and reaches 2.73 x 10^^ erg s“', compared with 
2.55 X 10^^ erg s ' in the baseline model. Also noticeable in 
is the decrease in the pre-breakout pulse associated with the 
deleptonization of the core on infall as described in §3. Shown 
here at f ^ 0.007 s, the pre-breakout peak reaches 0.75 x 10^3 
erg s“' with neutrino-electron scattering and only 0.52 x 10^3 
erg s“' without. By downscattering neutrinos during collapse 
and bounce via i'ee~ scattering, more neutrinos are able to es¬ 
cape before the core becomes opaque. Thus, there is a trade-off. 
Including inelastic scattering off electrons makes the collapse 
luminosity peak higher, but the main breakout luminosity peak 
lower. Although the breakout phenomenon is slightly affected, 
the spectra of Vg and Vg neutrinos at 200 ms are dominated by 
the charged-current interactions and so the effects of ve~ scat¬ 
tering on {Eu,), and (epj are modest, neutrinos are 

more significantly affected. Throughout the evolution, is 
decreased by ~30% by including v^e~ redistribution. is 
also profoundly altered. Comparing Fig. 21 with Fig. 7, we see 
that (ep^ ) is ~ 6 MeV higher without v^e~ redistribution. 

Figure 2 shows that there is significantly more deleptoniza¬ 
tion of the core at bounce with inelastic neutrino-electron scat¬ 
tering (thick solid line) than there is without it (long dashed 
line). The magnitude of the difference in core Yg at bounce 
is in keeping with other recent models (cf. Rampp 2000, §3). 
The fact that 7“"^® is higher at bounce in the latter case has im¬ 
portant consequences for the shock energy and the early post¬ 
bounce neutrino characteristics. In the model without neutrino- 
electron scattering, the peak positive velocity after bounce is 
higher (~2.5xl0^ cm s“', compared with ~2.0xl0® cm s“^) 
and peaks at larger interior mass (0.85 Mq, compared with 
0.7 Mq). In both models, as the shock stalls and the briefly pos¬ 
itive velocities achieved become negative, the matter behind the 
shock (the protoneutron star) oscillates on 5-10 millisecond 
timescales. As the protoneutron star pulsates, pressure waves 
move through the neutrinosphere of each species. These waves 
act to modulate the emergent luminosity and average neutrino 
energy. The amplitude of these near-periodic variations in neu¬ 
trino luminosity and energy are damped on ~10ms timescales. 
They are smaller for the model with neutrino-electron scattering 
because the shock in this model has less energy and the corre¬ 
sponding post-bounce oscillatory mass motions (‘ringing’) are 
thereby smaller at the neutrinosphere for each species. The 
smaller the local changes in temperature and density near the 
neutrinosphere, the smaller the changes in the local emission 
and absorption. Hence, we should expect smaller amplitude 
temporal modulations in the emergent neutrino spectral charac¬ 
teristics. 

Also potentially important is the fact that rises later in 
the model with neutrino-electron energy redistribution than in 
the model without such redistribution. Twenty milliseconds af¬ 
ter bounce. Fig. 20 shows that with scattering is ap¬ 
proximately half that for the model without this thermalization 
process. This has important implications for the detection of 
the Vg breakout pulse in terrestrial light-water Cerenkov neu¬ 


trino detectors like SuperKamiokande. In §8.1, we show that 
in such detectors, if one is to identify the Vg breakout signature 
uniquely, one must do so within 5-20 milliseconds of bounce, 
but before the Vg signal becomes appreciable. Hence, because 
neutrino-electron scattering suppresses the Vg signal during this 
epoch, it has consequences for the observation of the next galac¬ 
tic supernova. 


8. NEUTRINO DETECTOR SIGNATURES 


A total of just 19 neutrino events were observed from su¬ 
pernova SN1987A in the IMB (Irvine-Michigan-Brookhaven) 
and Kami oka II detectors. The signal lasted ~10 seconds and 
confirmed the basics of our theoretical understanding of core 
collapse; the amount of energy radiated in neutrinos was com¬ 
parable to that expected in protoneutron star cooling models at 
the time (Burrows & Lattimer 1986) and the inferred Vg neu¬ 
trino spectra had average energies of ^10-15 MeV (compare 
with Fig. 7). 

Far surpassing the sensitivity of Kamioka II and IMB, the 
current generation of neutrino detectors might collectively see 
thousands of neutrino events from the next galactic supernova 
(Burrows, Klein, & Gandhi 1992). Such an event would pro¬ 
vide exquisite neutrino spectral information and would put 
strong constraints on the physics of the core during collapse 
and explosion. Besides elucidating supernova physics, much 
information on the mass hierarchy of neutrinos could also be 
collected. For example, the neutrino signature might clearly re¬ 
veal whether a black hole or a neutron star had been created 
(Burrows 1988) and if the former is born, direct eV-scale mea¬ 
surements of v^- and i/ 7 .-neutrino masses could be made by the 
current generation of detectors (Beacom, Boyd, & Vogel 2000). 


Using the models presented in the preceding sections, we 
have made predictions about the dynamics and emergent neu¬ 
trino spectra expected from a class of supernova progenitors, 
with reasonable variations in microphysics and the models. In 
this section, we go a step further and fold our neutrino spectra 
with the thresholds and sensitivities of a subset of modern neu¬ 
trino detectors. The focus is on the detectability of the breakout 
phenomenon itself. 


In underground light-water detectors, such as a fully-repaired 
Super Kamiokande (SK), we can expect ^5000 neutrino events 
from the next galactic supernova at a distance of lOkpc. In such 
detectors the primary mechanism for neutrino detection is the 
Vgp ne^ process. The positron secondary emits Cerenkov ra¬ 
diation. For an electron or positron secondary, the differential 
number detected at distance D can be written (Burrows, Klein, 
& Gandhi 1992): 


dN = 


^ ^ da(e^,eg) 
A-ttD^ dSg 


X E{£g)S{t-t'- At)dt'dsgde^dt, 


(3) 


where Sg is the electron or positron energy, E(£g) is the detec¬ 
tion efficiency, Nj is the number of targets, L™™ is the number 
luminosity of neutrinos, f' is the source time, t is the detector 
time with D/c subtracted, Af(~ DI2c(mi,c^ j is the time 
delay of a neutrino with mass m^, and dajdSg is the differen¬ 
tial cross section for neutrino capture and electron or positron 
production. For most purposes dajdeg ~ a{e,y)5{e^-Eg- A), 
where A is the reaction threshold (in the case of the reaction 
Vgp ne^, A = m„-mp + mg ~ 1.804 MeV). Taking this ap- 
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proximation for da jclee, £(£«) = 1 , and = 0 , we obtain 


dN 

dt dSi, 


Nt 




(4) 


Using eqs. (3) and (4), it is a relatively simple matter to fold our 
time-dependent spectra with terrestrial detector characteristics, 
including efficiency and threshold corrections, in order to ob¬ 
tain the observed neutrino signal. We consider three different 
detectors, focusing on the ^g-rich breakout phase and the early 
neutrino signal that might be observed from the next galactic 
supernova. Although eq. (3) includes the effects of time-delay 
due to a finite neutrino mass, in what follows we assume mass¬ 
less neutrinos and employ eq. (4). 


8.1. SuperKamiokande 


The SuperKamiokande (SK) neutrino observatory is a light- 
water Cerenkov detector whose baseline volume for supernova 
neutrino detection is '^32 ktonne (Totani et al. 1998). We take 
as a threshold ~5 MeV (Beacom & Vogel 1998). The back¬ 
ground is typically of order several xO.l s“^ if the full detector 
volume is considered (Beacom & Vogel 1998) and will not fig¬ 
ure into our discussion here. 


Light-water Cerenkov detectors offer several different chan¬ 
nels for detection of the various neutrino species. The dominant 
reaction is charged-current absorption of Vg neutrinos on free 
protons (vep ne^). The positron secondary emits Cerenkov 
radiation, which is detected directly. Neutrino scattering on 
electrons, ^ also contributes for all neutrino species. 
For this process, we use the neutrino-electron scattering cross 
section (Sehgal 1974): 


dat 1 Uo 


dSg 2 MgC^ 

1 A 

cr,- = -^O-oAi 


A,+ ( 1-— ) Bi 


Sjy ^ 1 


nigC^ 2 


(5) 


where 


. 1 / 1 
Ai - - lAi+-Bi I , 


Ai = (Cv + CAf, Bi = (Cv-CAf 


and Cy = ± 1/2-t 2 sin^ 0 ^ 1 / for electron types and muon types, 
respectively. Ca = -l-l /2 for ly^ and and Ca = -1 /2 for Pg and 
Note that the neutrino-electron scattering cross section is 
6-7 times less for the /i- and r-type neutrinos, owing to the 
fact that for them the reaction can proceed via only the neutral 
current. Other important reactions include 


Vg+^^O- 

y^^¥+e- (eth = 15.4MeV), 

(6) 

Ug+^^O- 

y^^F+e- (eth = 1.66MeV), 

(7) 

Pe+'^O- 

>^^N + e^ (eth= 11.4MeV), 

(8) 

Vi+^^O 

^ly'i+y+X, 

(9) 


where e* is the energy threshold for the reaction. For the 
neutral-current excitation of (reaction 9), the 7 secondary 
has an energy in the range 5-10 MeV and, hence, is detectable 
by SK. Note that energy resolution for the secondary electron 
or positron in reactions ( 6 ), (7), and ( 8 ) is good to no better 
than ~ 10-20% (Burrows, Klein, & Gandhi 1992; Totani et 
al. 1998). For the study carried out here, we calculate event 


rates in SK due to only Pgp ne*' and neutrino-electron scatter¬ 
ing for all neutrino species. We assume 100% efficiency above 
the detector threshold. 

Figure 22 shows the integrated number of Vg and Pg neutrino 
events (logjQ[A(< t)]) detected in SK for our baseline IIMq 
progenitor (thick lines) and for the 20 Mq progenitor (thin 
lines) described in §5. The Vg neutrino signal from r'g-electron 
scattering (labeled, ‘vge~ Vge~') is the solid line for both pro¬ 
genitors. The Pg neutrino signal from charged-current absorp¬ 
tion on free protons (labeled, "Pgp ne"^’) is the dashed line 
for both progenitors. Note that many hundreds more Pg events 
accumulate from the 2 OM 0 progenitor (^ 1250) than from the 
11 Mq (~ 650) progenitor in the 250 ms of post-bounce evo¬ 
lution shown here. This important difference results from the 
steeper density gradient outside the iron core in the 11 M© pro¬ 
genitor and the correspondingly lower accretion luminosity (see 
Fig. 14). The inset in Fig. 22 shows the first 55 ms of the signal 
and the initial rise of the number Vg neutrinos detected due to 
breakout. The inset shows that approximately 5 Vg events can 
be expected from the collapse phase of the supernova, while 
10-15 Vg events might be observed solely from the breakout 
neutronization burst. This result seems to be independent of the 
mass of the progenitor. Although we do not include the results 
for the 15 Mq progenitor (see §5) in Fig. 22, as indicated by 
Fig. 13, its Vg breakout signature is virtually identical to that 
for the 11 Mq model. However, for the 15 Mq progenitor the 
luminosity of Pg neutrinos is slightly higher after 100 ms. As a 
result, 250 ms after bounce the integrated number of Pg events 
is ~5% higher than for the 11 Mq model. 

Note that the combined neutral-current and ^ signal 
due to scattering off of electrons is much smaller than that due 
to the process Vge~ Vge~. This signal never contributes sig¬ 
nificantly to the total neutrino event rate. In particular, less than 
1 event due to p- and r-neutrino scattering on electrons should 
be observed in the first 50 ms. 

Importantly for models of the neutrino oscillation phe¬ 
nomenon, the neutrino signal through collapse (~5 Vg events), 
breakout (10- 15 Vg events), and just post-breakout (^5 Vg 
events) should be composed of only Vg neutrinos. The break 
in the total number of events observed, which must occur as 
the Pg neutrinos begin to dominate the signal just 20-30 ms af¬ 
ter bounce, signals the end of this pure Vg neutrino phase. In 
addition, because the reaction Vge~ Vge~ is very forward- 
peaked and the process Pgp ne^ is slightly backward-peaked 
for £p^<15MeV (Vogel & Beacom 1999), one might hope to 
distinguish the early Pg signal from the post-breakout Vg neu¬ 
trino signal based on the directionality of the Cerenkov cones 
left by the positron and electron secondaries. However, if there 
is large mixing between Vg and or Vr neutrinos, the early Vg 
signal might be lessened significantly in SK (e.g. Takahashi et 
al. 2001 ). 


8.2. The Sudbury Neutrino Observatory 

The Sudbury Neutrino Observatory (SNO) in Sudbury, On¬ 
tario, is a light- and heavy-water neutrino detector. It is a spher¬ 
ical acrylic vessel 12 meters in diameter, containing 1 ktonne of 
D 2 O, that is surrounded by a cavity filled with light water. The 
fiducial light-water mass is approximately 1 .6 ktonne. SNO is 
most sensitive to and P^ neutrinos via breakup reactions on 
deuterons: 

Vi + d ^ n +p+Ui = 2.22M&N) ( 10 ) 
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Vi + d ^ n +p + Vi {e(a = 2.22Me,'W). ( 11 ) 

The liberated neutron secondary is thermalized in the heavy 
water and is then detected via capture on other nuclei within 
the detector volume. Currently, the neutrons from neutral- 
current events are detected via deuteron capture with an effi¬ 
ciency of ~30% (Ahmad et al. 2002; SNO Collaboration). Sev¬ 
eral schemes are currently proposed and in development for im¬ 
proved neutron-capture efficiency. The first is to lace the D 2 O 
with two tonnes of NaCl. The thermal neutron absorption cross 
section on ^^Cl is large (^83% efficiency) and would result in 
a 7 cascade peaked at 8 MeV. The 7 s would be detected by 
the ^9600 inward-looking photomultiplier tubes outside of the 
acrylic vessel. The second proposed neutron detection method 
involves hanging ^He proportional counters in the acrylic ves¬ 
sel and these would detect the neutrons directly (SNO collab¬ 
oration webpage^). Because the neutrons would have a '^5 ms 
mean capture time if the ^^Cl capture mechanism is employed 
(which would spread the turn-on), the neutral-current signal 
might be partially smeared with the ^^Cl approach. The neutral- 
current deuteron breakup reactions (reactions 10 and 11 ) also 
contribute to the Vg and Vg signal in SNO, but for them the 
charged-current deuteron capture processes, 

Ve + d ^ p + p + e~ {eik=lAAMeSf) ( 12 ) 

Pe-tc/^ n-l-n-l-e'^ (eth = 4.03MeV), (13) 

are more important above e^, ~ 10 MeV (Burrows, Klein, & 
Gandhi 1992). Process (13) can be used to observe Vg neutrinos 
uniquely via the simultaneous detection of the secondary neu¬ 
tron and the Cerenkov emission from the final-state positron. 
Although we include the threshold energy in our calculations 
of the total signal, as with SK we assume a 100% efficiency for 
SNO above the detector theshold. 

In Fig. 23 we show the integrated number of Vg neutrinos 
detected in SNO via the processes Vgd ppe~ (dotted line), 
Vgd npvg (solid line), and Vge~ Vge~ (dashed line) in the 
first 250 ms after bounce in our IIMq baseline model. The 
total Vg signal is the thick solid line. For comparison, we pro¬ 
vide the contribution to the total number of detected events 
from all other neutrino species (long dashed line, labeled ‘To- 
tal(other)’). This signal includes the Vgd npVg, Vgd nne^, 
and Vgp ne^ processes in the light-water portion of the detec¬ 
tor, as well as the processes v^.rd npv^ T, npv^ T, 

and neutrino-electron scattering throughout the entire detector 
volume. Note that it might be possible to exclude the charged- 
current absorption process Vgd —> nne*' by cutting on events 
that have simultaneous neutron and positron detection. This 
would lower somewhat the line labeled Total(other) and make 
the early i^g breakout signal easier to see. However, such a pro¬ 
cedure requires accurate signal timing, which would be a chal¬ 
lenge if the ^^Cl neutron capture mechanism is employed. As it 
stands, one might expect to get ~5 events in SNO from the i^g 
breakout burst in the first ~10 ms after bounce. 

8.3. ICARUS 

Although it has not yet reached its goal mass, ICARUS"^ is 
designed to be a 3.6-ktonne drift chamber of pure liquid "^®Ar 
that will be sensitive primarily to electron neutrinos through 
electron-neutrino capture on Ar: 

^ http://www.sno.phy.queensu.ca/sno/sno2.html#nc 
^ http://www.aquila.infn.it/icaras/ 


Recent shell-model calculations have shown that this process 
can proceed through super-allowed Fermi transitions to the 
4.38 MeV excited isobaric analog state, as well as through 
several Gamow-Teller (GT) transitions to other lower-lying 
states in (Ormand et al. 1995). This has recently been 
investigated experimentally by studying the /^'''-decay from 
‘^‘’Ti (Bhattacharya et al. 1998). Employing isospin symmetry, 
Bhattacharya et al. (1998) inferred the transition strengths of 
“^^Ar K from '^^Ti Sc. These theoretical and experi¬ 
mental results indicate that the GT transitions enhance the total 
i^g capture cross section by a factor of three over the pure Fermi 
transition cross section. We take the super-allowed Fermi tran¬ 
sition cross section from Raghavan (1986); 

ctf-1.02 X 5.365)^ cm^ (e* = 5.885MeV) (14) 

and we assume that the total lyg absorption cross section is given 
by = 3crF. 

The ICARUS detector is currently at 600 tonnes and has re¬ 
cently been approved to expand to 3 ktonnes. Because of its 
sensitivity to I'g neutrinos, ICARUS is of particular importance 
in detecting the early breakout pulse so prominent in Fig. 6 . In 
the results presented below, we assume that the detector effi¬ 
ciency is 100 % above threshold and that the detector mass is 
3 ktonnes. 

The I'g breakout signal in ICARUS is very likely detectable. 
Figure 24 shows N(< t) for the ICARUS detector (with 3 
ktonnes of Ar). We plot the two dominant detection chan¬ 
nels for Vg neutrinos: VgAr absorption (thick solid line) and 
Vie~ Vie~ (thin solid line). We also include the event rate, 
dN/dt (in units of 50 s“^), so that the width of the Vg breakout 
spike can be compared with the detector signal in the very early 
phase. Note that the other neutrino species will contribute to the 
total signal through only the Vie~ process and that the 

combined signal amounts to only two events in the first 250 ms 
shown in Fig. 24. Combining the r'^-electron scattering and Vg- 
Ar absorption rates, ICARUS should expect to detect approxi¬ 
mately 10 Vg events from the breakout pulse alone, more than 
SNO and as many as Super K. 


8.4. The Influence on the Signal of Various Neutrino 
Processes and the EOS 

8.4.1. Weak Magnetism & Recoil 

The weak-magnetism/recoil correction is relevant both in the 
supernova itself and in the detection of Vg neutrinos in light- 
water detectors via the process VgP ne'^. In §7.2 we com¬ 
pared the resulting spectra from models of core-collapse with 
and without this correction. In all our detection calculations we 
retain the weak-magnetism/recoil correction in calculating the 
signal due to Vgp ne^. Comparing the detector signal with 
and without this correction we find differences in events rates 
for neutrinos of all flavors of less than 1%. The models that 
include the correction have slightly higher event rates. For anti¬ 
neutrinos, in particular for Vg in SK, we find that the event rate 
with the correction is larger by approximately 7%. 
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8.4.2. Inelastic Neutrino-Electron Scattering 

Inelastic neutrino-electron scattering in our core-collapse 
models significantly affects the observed neutrino signal. 
The total event rate in SK via neutrino-electron scatter¬ 
ing 200 ms after bounce is ^15% higher in the model with¬ 
out inelastic neutrino-electron scattering in the supernova. If 
this process is neglected in the supernova models, the signal 
in SNO due to neutral-current deuteron breakup by neutri¬ 
nos is nearly doubled. Approximately 200 ms after bounce, the 
corresponding event rates are 110 s“' and 65 s“', respectively. 

Inelastic neutrino-electron scattering also affects the ob¬ 
served early-time D,, signal in SK. Figure 20 demonstrates that 
the Dg luminosity very near breakout is larger in the model with¬ 
out neutrino-electron scattering, though at later times models 
with and without neutrino-electron scattering are quite similar. 
At 200 ms after bounce, we obtain in SK 50 more Dg events in 
the model without neutrino-electron scattering (550 versus 500 
events). Just 20-30 ms after bounce, the event rates differ by 
10-15%. 


8.4.3. Nucleon-Nucleon Bremsstrahlung 

Including nucleon-nucleon bremsstrahlung as a production 
process for neutrinos softens and brightens the resulting 
spectrum (see §7.3 and Fig. 18). By including bremsstrahlung, 
the total event rate in SNO due to neutral-current deuteron 
breakup by neutrinos 200 ms after bounce is increased by 
13%. Similarly, the total number of neutrino-electron scatter¬ 
ing events in SK is increased by 20% (10.5 versus 8.5 events 
after 200 ms). The event rates during this epoch of accretion- 
driven neutrino luminosity are 62 s“^ and 55 s“', respectively - 
a fractional difference of ^ 11 %. 


8.4.4. LSEOS 

In § 6 , we presented results from core-collapse calculations 
employing different nuclear compressibilities in the context of 
the LSEOS (see Appendix A). In the model with k = 180MeV 
we find that the r'^-electron scattering event rate in SK 200 ms 
after bounce is 160 s“' and that the corresponding cumulative 
event total is 47. The model with k = 375 MeV gives 150 s“^ 
and 42, respectively. The event rates in SK due to Vg capture on 
free protons 200 ms after bounce are 2870 s“' and 3080 s“' for 
K = 375 MeV and k = 180 MeV, respectively - a fractional dif¬ 
ference of just ~ 7%. The signal in both SK (due to neutrino- 
electron scattering) and in SNO (due to deuteron breakup) are 
modified at the 10 % level, with the larger event rates corre¬ 
sponding to the lower nuclear compressibility (k = 180MeV). 


8.5. Neutrino Oscillations 

Due to severe matter suppression of standard flavor neu¬ 
trino oscillations at the high densities encountered in supernova 
cores, such oscillations do not appreciably affect the dynamics, 
neutrino spectra, nor neutrino luminosities in the interior. We 
have calculated a few models that approximately incorporate 
flavor oscillations employing the oscillation parameters implied 
by the solar neutrino and atmospheric neutrino deficits and find 
almost no effect on core evolution. 


The standard Boltzmann equation involves intensities, not 
amplitudes and phases. It is the evolution equations for the lat¬ 
ter, or for a full density matrix, that should be calculated in 
order to take neutrino oscillations properly into account. How¬ 
ever, we have opted for a simpler approach that captures the 
basic effect. For a freely propagating electron neutrino, the 
probability that it oscillates into a neutrino is: 



where 9m and Em are the oscillation angle and length in matter, 
r is the distance propagated, and 9 is the oscillation angle in 
vacuum. Em and 9m are given by the expressions: 

• 2^n sm^29 

Sin 28tr 


sin^ 29-\-(cos29-Ly/L q)^ ’ 


and 

Em “ Ey 
where 


1-?^Ecos20 + (^)2 
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1 - 1/2 


= E„ 


sin 20m 

sin 20 


E,= 


4TTJiei, 

drn^c 
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V 10 MeV/ 


5x10-" eV 
5m^ 


(16) 


(17) 


(18) 


is the vacuum oscillation length, Lq ~ 3.3 x cm/pn, 

P \2 = p/lO'^g cm"^, and 6rn^ is the standard mass-squared- 
difference. The matter suppression of neutrino flavor oscil¬ 
lations is quite severe. At resonance, 9m = 7 r /4 and E,,/Lo = 
cos 20 , but for the range of densities that obtain in protoneu¬ 
tron stars and the values of 5m^ (~ 5 x 10 -"eV^) and sin ^20 
(> 0.8) derived from the KamLAND (Eguchi et al. 2003) and 
solar neutrino data, there are no MSW resonances in supernova 
interiors. We have not taken into consideration the possibil¬ 
ity of density fluctuations induced neutrino flavor depolariza¬ 
tion (Loreti et al. 1995). If we assume that a hard scattering 
or absorption decoheres the state (“resets the clock”), then a 
rough measure of the probability of oscillation from Vg to in 
a radius interval dr is Pe^e~''l^dr/K, where A is the extinction 
mean-free path. The appropriate integral gives us the approxi¬ 
mate transition rate 

rg^ =——sin^ 9m (—) ■ (19) 

There are corresponding expressions for F^^, etc. These rates 
are plugged into rate equations that conserve total neutrino 
number: 


^-^ = Tg^J,-T^gJ,jA, (20) 

where 7, is the energy dependent zeroth moment of the specific 
intensity for neutrino species / (x 4 for the v^A). Hence, using 
eqs. ( 20 ), we can incorporate the effects of neutrino oscillations 
on the internal supernova dynamics and spectra. As stated, they 
are small. 


The edge of our computational grid resides deep inside the 
supernova progenitor. The spectra we have calculated and 
folded with terrestrial neutrino detectors will undoubtedly be 
modified by oscillation effects in propagating from the edge 
of our computational domain to the detector through the low- 
density progenitor envelope. Hence, our neutrino number spec¬ 
tra must be multiplied by effective survival probabilities tak¬ 
ing into account the density profile of the overlying progenitor 
and, depending upon the position of the detector, the strata of 
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the Earth itself. Because our focus is here on producing re¬ 
liable emergent spectra from core-collapse and the diagnosis 
of the basic effects on this signal of changes in the supernova 
microphysics, and because modifications to the signal in detec¬ 
tors depends on the wide range of progenitor envelope profiles, 
oscillation parameters, and Earth day-night effects, we save a 
detailed investigation of neutrino oscillations for a future work. 
The reader is referred to Ando & Sato (2002), Dighe (& Smirnov 
(2000), Dutta et al. (1999), and Takahashi et al. (2001) for re¬ 
cent discussions of possible neutrino oscillation effects in the 
detected supernova neutrino signal. 

9. SUMMARY AND DISCUSSION 

We have constructed fully dynamical models of core col¬ 
lapse in spherical symmetry, focusing on the shock breakout 
phenomenon and the first lOO’s of milliseconds after bounce. 
We employed a newly-developed algorithm for radiation- 
hydrodynamics, which gives a Boltzmann solution to the neu¬ 
trino transport problem using the Eeautrier technique, the 
tangent-ray method, and Accelerated Lambda Iteration. The 
code provides good angular and spectral resolution of the ra¬ 
diation field and incorporates realistic neutrino microphysics. 
We find that at least 20 energy groups are required to ade¬ 
quately resolve the neutrino spectrum and the dynamics; 10 en¬ 
ergy groups are not enough. Eurthermore, the tangent-ray algo¬ 
rithm provides very good and adaptive angular resolution that 
is well-suited to the spherically-symmetric supernova problem 
in which the neutrino fields transition from a diffusive interior 
to a free-streaming exterior that nevertheless still resides within 
the relevant computational domain. Using this scheme, we have 
presented some of the first high-resolution angular distribution 
functions for neutrinos in stellar collapse. 

We have explored the effects on dynamical simulations of 
core collapse and its neutrino signatures of nucleon-nucleon 
bremsstrahlung, inelastic neutrino-electron scattering, weak- 
magnetism/recoil, and artificial decreases in the neutrino- 
nucleon scattering rates. Our new algorithm for incorporating 
inelastic scattering explicitly as a source and a sink in the colli¬ 
sion term of the Boltzmann equation is both efficient and robust 
(§B). Eurthermore, we have studied the dependence of the col¬ 
lapse dynamics and neutrino signal on variations in the equation 
of state, as parametrized by Lattimer and Swesty (1991). Sur¬ 
prisingly, the dependence of the salient quantities of collapse 
on the EOS (again, as parametrized by Lattimer and Swesty) is 
weak. 

We have computed the collapse, bounce, and post-bounce 
evolution and neutrino signatures for three progenitors: 11, 15, 
and 20 Mq. In the process, we have compared key quanti¬ 
ties such as core entropy and electron fraction at bounce, peak 


neutrino breakout luminosity, and the evolution of the neutrino 
spectra over hundreds of milliseconds with the corresponding 
quantities in other recent supernova simulations that employ 
comparably sophisticated algorithms. We evolve most models 
to several hundred milliseconds after bounce and although all 
of our models form regions of net neutrino heating behind the 
shock, none yields an explosion in the first 250 milliseconds. 

We have folded our neutrino spectra with the sensitivities and 
thresholds of several underground neutrino detectors and have 
focused on those detectors most likely to observe and uniquely 
identify the electron-neutrino breakout burst. We show that 
SuperKamiokande, the Sudbury Neutrino Observatory, and the 
ICARUS detector might all see this important signature of core 
collapse. We find that the collapse and breakout electron- 
neutrino signal in these detectors is remarkably similar (see 
Eig. 22) for all three progenitors considered in this work. The 
pre-collapse profiles are different in entropy and composition, 
but the breakout peak is roughly the same in each model. These 
results are further confirmed in preliminary runs with a TOM© 
progenitor. 

Eigure 22 shows that the event rate due to charged-current 
capture of neutrinos on free protons for the 20 M© progen¬ 
itor (5850 s“') is more than double that for the 11 M© model 
(2920 s“') 175 ms after bounce. The nearly linear rise \nN{<t) 
reflects the fact that the accretion luminosity is roughly con¬ 
stant in this early post-bounce epoch. The difference between 
the 11 and 20 M© models can be attributed primarily to a dif¬ 
ference in M (see Eig. 14), which, in turn, is a consequence of 
the relatively steeper envelope density profile of the lighter pro¬ 
genitor. The behaviors of the density profiles of the 11 M© and 
20 M© models roughly bracket those found in the full progeni¬ 
tor model suite of Woosley & Weaver (1995). Hence, we expect 
that most other progenitors should yield accretion luminosities 
that produce signals between those of the 11 M© and 20 M© 
models in Eig. 22. 
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APPENDIX 

TABULAR NUCLEAR EQUATION OF STATE 

We employ the equation of state due to Lattimer & Swesty (1991) (the LSEOS), based on the finite-temperature liquid drop model 
of nuclei developed in Lattimer et al. (1985). We have constructed an efficient three-dimensional tabular version of this EOS, taking 
the matter temperature (T), mass density (p), and electron fraction (Ye) as variables. The table consists of 180 equally spaced points 
in log[Q[p (g cm“^)] with 6.4 < log[Q[p (g cm“^)] < 15.1, 50 equally spaced Ye planes with 0.05 <Ye < 0.51, and 180 equally spaced 
zones in log[Q[7’ (MeV)]. The table is not cubic. Eor log[Q[p (g cm“^)] = 6.4, -0.8 < log[Q[7’ (MeV)] < 1. At log[Q[p (g cm“^)] = 15.1, 
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0 < log[Q[r (MeV)] < 1.6. We find this range of temperatures and densities sufficient for all calculations we perform. At each point 
in the table we save the specific internal energy (E), the pressure (P), the entropy per baryon (s), the specific heat at constant volume 
(Cv), the neutron, proton, alpha particle, and heavy nucleus fractions (A„, Xp,Xa, Xh, respectively), A and Z/A for the heavy nucleus, 
the electron chemical potential (/Xg), = dlnP/d\np\s, dP/dT\p, /I = pn — Pp, and the derivatives dp/dT and dp/dYg. In total, our 
tabular LSEOS takes up ~200 Megabytes of memory. Given T, p, and Yg, the EOS performs three six-point bivariant interpolations 
(Abromowitz & Stegun 1972) in the T -p planes nearest to and bracketing the given Yg point. A quadratic interpolation is then 
executed between Yg points to obtain the desired thermodynamic quantity. This procedure is employed for all quantities in the table 
except the mass fractions of neutrons, protons, alpha particles, and the heavy nucleus. Eor these quantities, two four-point bivariant 
interpolations in the T -p plane are combined with a linear interpolation between Yg planes. The table uses integer arithmetic to 
find nearest neighbor points, thus alleviating the need for time-intensive search algorithms. Because most hydrodynamics routines 
update specific internal energy, we include a Newton-Raphson/bisection scheme which iterates on temperature, given a fixed internal 
energy, until the root is found to within a part in 10*. Similar iteration routines are employed if one wishes to iterate on entropy or 
pressure. 

The LSEOS extends down to only ~ 5 x 10® g cm“* and its validity in this density regime is guaranteed only for fairly high 
temperatures - the assumption of NSE being thereby satisfied. Eor calculations in which an explosion occurs, the shock will quickly 
evolve down the progenitor density gradient to regions where the LSEOS breaks down. Eor this reason we have coupled to our tabular 
version of the LSEOS the Helmholtz EOS (Timmes & Arnett 1999; Timmes & Swesty 2000), which contains electrons and positrons 
at arbitrary degeneracy and relativity, photons, nuclei and nucleons as non-relativistic ideal gases, and Coulomb corrections. At p = 
6 X 10^ g cm“* we assume the LSEOS is valid. At p = 4 x 10^ g cm“* we employ only the Helmholtz EOS. Eor densities between the 
two, we quadratically interpolate all relevant thermodynamic quantities. Because the LSEOS assumes NSE and the Helmholtz EOS 
takes non-NSE abundances, there is no thermodynamic consistency between the two EOSs. Eortunately, in thermodynamical regimes 
relevant for core collapse at these densities, the electron-positron/photon component of the matter dominates the pressure. Therefore, 
one expects few hydrodynamic artifacts in piecing together these two equations of state. In fact, given a realistic composition for the 
Helmholtz EOS, at a density of p = 5 x 10^ g cm“*, pressure differences between the LSEOS and Helmholtz EOS are ^ 1% in the 
stellar profiles we employ in this study. Differences in entropy, however, are of order 5-10%. 

We have performed an extensive set of tests of our tabular implementation of the LSEOS to ensure that thermodynamical con¬ 
sistency is maintained during dynamical simulations. These tests include the adiabatic compression of a single fluid element over 
eight orders of magnitude (from an initial density of 10^ g cm“*) and the purely hydrodynamical core collapse of realistic supernova 
progenitors. Eor the single-element tests, we use the fractional change in the element entropy (As/si, where Si is the initial entropy 
at low density) to assess the consistency of the LSEOS as the element is compressed and the density increases to ~ 5 X 10'“* g cm *. 
Eor fractional changes in the density (or specific volume) at each step of 1%, our table gives As/si ~ 0.5%. This is to be compared 
with the results using the actual analytic LSEOS, from which the table was composed, which yield As/si ~ 0.25%. The number 
of temperature and density points used in the table does effect Ai/si. Eor example, using just 100 T, 100 p, and 50 Yg points over 
the same ranges as the larger table raises As/si to about ~ 1.25%. In all our calculations we have opted for the larger table. Eor 
the adiabatic collapse tests, we coupled the LSEOS to a variety of one-dimensional hydrodynamics schemes. With radiation turned 
off and before any shocks form, each mass zone should remain isentropic during the collapse of a massive progenitor star. Using 
a variety of zoning schemes with both artificial-viscosity and Riemann solvers for shock capturing and resolution, we found that 
higher-order EOS interpolation schemes as described in Swesty (1996) are not needed to maintain the degree of adiabaticity (1%) 
during collapse and bounce that we were able to achieve with our more straightforward tabular and interpolation routines. 


INELASTIC NEUTRINO SCATTERING 


A number of authors have addressed the issue of inelastic neutrino-electron scattering as an important energy redistribution process, 
which helps thermalize the neutrinos and increases their energetic coupling to the matter in supernova explosions (Bruenn 1985; 
Mezzacappa & Bruenn 1993). Here, we review the Legendre expansion formalism for approximating the angular dependence of the 
scattering kernel, and detail our own implementation of scattering terms in the Boltzmann equation. 


The general collision integral for inelastic scattering may be written (Bruenn 1985) as 




c(27r7ic)* 

scatt 


/'R”(e^,e'^,cos6»)-/i, 


c(27r7ic)* 




t _ -scatt r 

Ku Ji 


(Bl) 


where cos 0 is the cosine of the scattering angle, e,g is the incident neutrino energy, and is the scattered neutrino energy. Although 
we suppress it here, the incident and scattered neutrino phase space distribution functions (f,g and /^, respectively) have the following 
dependencies; f,g = f^ipt , p^e^) and = fl(r,t,p'p and p’ are the cosines of the angular coordinate of the zenith angle in 
spherical symmetry and are related to cos 9 through 

cos9 = pp + [(l- p\l- cos((j3 - (j)'). (B2) 

R” is the scattering kernel for scattering into the bin (e,^, p) from any bin (e[,, p') and is the scattering kernel for scattering 
out of the bin (e^, p) to any bin (e^, p'). The kernels are Green functions, correlation functions that connect points in energy 
and momentum space. In fact, one may also write R(£^,£l,cos9) as R{q,uf), where u;(= — is the energy transfer and q{= 

[£^-l-£j,^-2ei^£j,cos0]^/^) is the momentum transfer, so that the kernel explicitly reflects these dependencies. 



15 


It is clear from eq. (Bl) that the source, or the net scattering into the beam (e^, /i) is a function of the occupancy (//,) in all other 
beams. The sink term also depends on //, through the blocking term (1 -fl), reflecting the Fermi-Dirac statistics of the neutrino. 
A solution to the Boltzmann equation yields fi, (and, hence, //,) at all times, radii, energies, and angles, fi, is not known a priori 
and cannot be assumed to be Fermi-Dirac. Only in equilibrium should approach a Fermi-Dirac distribution, characterized by the 
local temperature and with a chemical potential that reflects the local neutrino number density. The transport problem is difficult 
enough without the added complication of non-linear integral source terms. The full implicit solution, including energy and angular 
redistribution, is numerically cumbersome. Instead, we make several simplifications that make the problem tractable, efficiently 
solved, and explicit. 

An important simplification comes from detailed balance, a consequence of the fact that these scattering rates must drive the 
distribution to equilibrium; R'" = where /3 = l/T. We may therefore deal only with R°'^\ The scattering kernels generally 

have complicated dependencies on scattering angle and the angular integration over scattered neutrino phase space, implicit in 
eq. (Bl), cannot be computed analytically. For this reason, we approximate the angular dependence of the scattering kernel with a 
truncated Legendre series (Bruenn 1985). That is, we take 


, cos 0) = ^ , cos 6»)P/(cos Q ), 


/=o 


where 


^i{£^,£'^) = J d(cos9)R°'"\£„,£l,cos9)Pi(cos9). 


In practice, we expand only to first order so that 


,cos9) ^<^o(£^,£'J+^<i’i(£^,£'J cos 9. 


Substituting into the first term on the right-hand-side of eq. (Bl) (the source) gives 
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Substituting for cos 9 using eq. (B2) and using the definitions 
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we have that 
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Integrating over /x to get the source for the zeroth-moment equation. 
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Similarly, employing the Legendre expansion, we may write the sink term of the Boltzmann collision integral in terms of the angular 
moments of /': 
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The contribution to the zeroth-moment equation is then 
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Combining these equations, we find that 
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One can see immediately that including another term in the Legendre expansion, i.e. taking 

f?™' ~ (1 /2)T>o -I- (3/2)T> 1 cos 61 -I- (5/2)T>2( 1 /2)(3 cos^ 6» - 1), 

necessitates including P^, and P'^, the second angular moment of the neutrino phase-space distribution function, in the source and sink 
terms. Our transport scheme porvides both the zeroth- and first-moment equations with the spectral Eddington factors p^, = Pu/Jv 
and gi, = Nu/Ju, where 


~ 1 

^-2 


c+\ 


1 


r+1 


dp,p, U and N^ = -f dpp^f^. 


(B13) 


It is therefore straightforward for us to include terms up to R°'^^ oc cos^ 9 in our Legendre expansion of the scattering kernel. In 
practice, however, we retain only the linear term. 
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Neutrino-Electron Scattering 

The scattering kernel R(e,y,el,cosO) in §B is related to the fully relativistic structure function for neutrino-electron scattering; 

9) = 2G^ [A5i {q,uj) + Siiq, to) + BSsiq, w)](l - . (B14) 

The relativistic structure functions (5,) are given in terms of the imaginary part of the polarization functions (If). Each of the retarded 
polarization functions can be written in terms of one-dimensional integrals over electron energy (e,.) (see Reddy, Prakash, & Lattimer 
1998). 

Even though we have made the simplifying assumption that the scattering kernel can be approximated by a Legendre series 
truncated at first order (eq. B5), in a fully dynamical calculation it is numerically costly to compute the Legendre moments of 
the scattering kernel (^o(£y,£'y) and ^i(ey,£l)) via eq. (B4) at each point on the radial grid (at each temperature, density, and 
composition point) at each time step. In order to decrease the computation time, we tabulate $o(ei/,eJ^) and <i)i(£^,£j,) for each 
and for each neutrino species, on a grid in temperature and At each T point, we calculate both kernels for each 
point, given the energy grouping for the particular calculation. The angular integrals over cos 9 for I = 0 and 1 = 1 in eq. (B4) are 
carried out using 16-point Gauss-Legendre quadratures. During an actual supernova calculation, we use simple trapezoidal rule 
quadrature to calculate the integral over e'^ for a given £y. Each term in this integral contains the kernels, which must be evaluated 
for each ey-e'y pair at the temperature/density/composition point currently being addressed, qe is calculated by the equation of state 
and we do a six-point bivariant interpolation in T-q^ space, for the given Ey-e'y combination. In practice, we use 40 energy groups 
(nf), 30 temperature points (Nt), and 30 qe points (Njj). The tables are then I x nf x n'jr x Nt Njj in size, with I = 2 ($o and <i)i), 
or approximately 50 Megabytes. Since the vector and axial-vector couplings for electron scattering are different for each neutrino 
species, we compute tables for lOe, Vg, and for each dynamical calculation. and are then included explicitly as a source 
and sink, respectively, in a manner analogous to any of the absorption and emission processes detailed in Burrows (2001). Using this 
method, our calculations including neutrino-electron scattering, are just 10-15% slower than our calculations ignoring this important 
equilibration process. 

Because we do not evaluate the scattering source and sink implicitly, we introduce an explicit timescale into the energy and electron 
fraction updates returned by our transport algorithm. In effect, if the scattering timescale (cx^"^^") * is shorter than our timestep, we 
may encounter a numerical instability. Eor this reason, at high densities, where the neutrino phase-space distribution function is in 
local thermodynamic equilibrium, we divide the source and sink by a factor (typically 100 above p = 10^"^ g cm“^). This reduces 
the rate artificially, thus increasing the timescale for scattering. As this process is totally negligible in this regime, particularly 
considering the fact that fy = ft^ at these high densities, this approximation is acceptable. 
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M (MJ M (M^) 


Fig. 1.— Velocity v (in cm upper left panel), log 2 Q[p] (in g cm“^, upper right panel), entropy s (in kg baryon”^, lower left panel), and temperature T (in MeV, 
lower right panel) in the IIMq progenitor, as a function of mass coordinate in Mq for five snapshots in time. The thin solid line shows the initial configuration 
and the thick solid line shows the model at bounce. The dotted, dashed, and dot-dashed lines are snapshots at 1.5 ms, 17 ms, and 200 ms after bounce, respectively. 
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Fig. 2.— The electron fraction Ye in the IIMq progenitor, as a function of mass coordinate in Mq for five snapshots in time. The thin solid line shows the 
initial configuration and the thick solid line shows the model at bounce. The dotted, short-dashed, and dot-dashed lines are snapshots at 1 ms, 17 ms, and 200 ms 
after bounce, respectively. For comparison, we include the bounce Ye profile with inelastic neutrino-electron redistribution turned off (long dashed line). Compare 
this figure with Fig. 1, which shows the basic hydrodynamical evolution for this same model. 
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Fig. 3.— Snapshots of in 10^"^ erg s“^ as a function of mass for the fiducial 11 Mq progenitor, showing the whole envelope of luminosities realized in this 

progenitor from bounce through electron-neutrino breakout. Times re relative to hydrodynamical bounce. The initial rise in at 0.3 Mq comes from electron 
capture on newly liberated free protons before the shock forms, but just after bounce. The second, larger peak forms after shock formation. The temperature and 
density increase across the shock dissociates nuclei into free nucleons. Subsequent electron capture on free protons generates the breakout pulse, which reaches 
Lt/g 1.5 X 10^“^ erg s“^ locally. 
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Fig. 4.— Temperature (MeV, dashed line) and electron neutrino luminosity (10^^/15 erg s“^, thick solid line) as a function of enclosed mass 3.2 ms after bounce. 
The four small insets (thin solid lines on dotted axes) are polar plots of the angular distribution of the electron-neutrino specific intensity (i.e., (0)) for = 12 MeV 

at four mass points: (A) 0.674 Mq, (B) 0.89 Mq, (C) 0.945 Mq, and (D) 1.03 Mq, each point indicated specifically by large dots on the temperature profile. At 
(A), the radiation field is isotropic, indicating that the flux is zero and that the neutrinos are trapped. At (D), the radiation field is nearly decoupled from the matter, 
as indicated by the forward-peaked angular profile. 
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Fig. 6.— Lue (thick solid line), (thin solid line), and (dotted line) measured at the outer edge of the grid in erg s“^ as a function of time for the fiducial 
M= 11 Mq progenitor. Time is measured relative to bounce. Note that we define r = 0 as the time of hydrodynamical bounce. The finite light travel-time to the 
edge of the grid creates a 7 ms offset between hydrodynamical bounce and the initial dip before the large Ug breakout pulse. 
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Fig. 8. — Luminosity spectrum of Ve neutrinos at infinity at various pre- (thin solid lines) and post-breakout (dashed lines) times. In this figure, time is measured 
relative to the peak breakout spectrum (thick solid line). The thin solid lines correspond to 11.6, 5.1, and 1.3 ms before the peak and the dashed lines denote the Ug 
spectrum 4.2, 9.6, and 40.5 ms after the peak. 
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Fig. 10.— Polar plot of the normalized specific intensity, constructed from the Feautrier variables at the neutrinosphere (as defined in eq. 1) for Su ~ 10.7 MeV 
for Ve neutrinos. This plot shows only every other energy grid point (solid lines). There are 269 angular bins in each quadrant. At the largest energies, the radiation 
field is nearly isotropic and the flux is quite small. At the lowest energies, the radiation field is beginning to decouple from the matter and 1^, is forward-peaked. 
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Fig. 11.— Polar plot of the specific intensity as a function of 0 for Ve neutrinos with Si/ = 10.7 MeV at selected radii. We maintain the same scale as in 
Fig. 10 for the line labeled ‘60 km’ and normalize the absolute value for all other radii to that reference 1^,. Thick solid lines denote lugiO^Eu = lO.TMeV) at 
the neutrinosphere ~ 60 km), at the gain radius {Rg ~ 100 km), and at the shock radius {Rs ~ 200 km) 200 ms after bounce. Thin solid lines show the specific 
intensity at various intermediate radii. 





Fig. 12.— The neutrinospheres (as defined by eq. 1) as a function of neutrino 
line) for four post-bounce times in the baseline 11 Mq model. 



for each neutrino species {Ue, thick solid line; Ug, thin solid line; dotted 
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Fig. 14.— Mass flux M in Mq s“^ as a function of radius for our 11 Mq (thick solid line), 15 Mq (thin solid line), and 20 Mq (dotted line) models, all at 
200ms after bounce. Open triangles mark the Ve neutrinosphere at {sue) — 12, 12, and 13.5MeV for the IIMq, 15M0, and 2OM0 models, respectively. 
Filled squares and open circles mark 2::: 15, 15, 16MeV) and ) — 20, 20, 21 MeV) for the 11 M0,15 M0, and 20 M0 models, respectively. 
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Fig. 15.— Li/g and in erg s ^ at infinity as a function of time, for the fiducial model and for the model with artificially decreased neutral-current neutrino 
neutron and neutrino-proton cross sections. 
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Fig. 18.— as a function of radius for 11 Mq models with bremsstrahlung and inelastic neutrino-electron scattering (solid line), without bremsstrahlung 

and with inelastic neutrino-electron scattering (short dashed line), and with bremsstrahlung and without inelastic neutrino-electron scattering (long dashed line), 
approximately 220ms after bounce. The inset shows the neutrinosphere Ru^, dehned by eq. (1) as a function of energy in the same models. 
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Fig. 19.— and Li/^ of the emergent spectrum at infinity as a function of time, for IIMq models with (thick solid line) and without (thin solid line) 

nucleon-nucleon bremsstrahlung. 
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Fig. 21.— {^ue) (thick solid line), {sue) (thin solid line), and {si^^) (dotted line) at infinity for the 11 Mq progenitor without inelastic neutrino-electron scattering. 
Compare with Fig. 20, which shows the corresponding luminosities and Fig. 7, which shows the average energy evolution in our fiducial model (including inelastic 
neutrino-electron scattering). 
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Fig. 22.— The integrated number of neutrino detection events (logjQ[A^(< r)]) in SK for our 11 Mq (thick lines) and 20 Mq (thin lines) models via electron 
scattering {i'ee~ — Uee~, solid lines) and i?e capture on free protons {Pep dashed lines) for supemovae at lOkpc. The inset is an expanded view of the 

i/e breakout signal. Note that both models imply that SK should see a distinct and observable Ue signature in the first 30 ms. Nearly 20 Ue events accumulate 
before being swamped by the dominant Pe signal. For clarity, we do not include the corresponding lines for the 15 Mq progenitor in this figure. The Ue neutrino 
signal in the 15 Mq model is identical to that of the 1 1 Mq in the first 50 ms. The Pe event rate is slightly larger than that for the 11 Mq model; the 15 Mq model 
accumulates 680 events whereas the 1 1 Mq model yields 650 events in the 250 ms after bounce shown here. 
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Fig. 23.— Integrated number of Ve neutrinos detected in SNO for a supernova at lOkpc as a function of time. The individual contributions to the detected signal 
via the Vgd —s- ppe~ (dotted line), Ved —s- npVe (solid line), and Vee~ —^ Vee~ (dashed line) processes are shown. The sum is the thick solid line. The sum of all other 
processes, including Ved —^ npVe, Oed nne^, Pep in the light-water portion of the detector, as well as P^.rd —> upp^^r, i^iJ.,Td ^ npu^^r, and Pe-, 

and F^,r-electron scattering throughout the entire detector volume, is shown as the long dashed line. This signal corresponds to our fiducial 11 Mq core-collapse 
model. 
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Fig. 24.— Integrated number of Ve events as a function of time via absorption on Argon (labeled V^Ar —s- K* +^“’, solid line), including the Fermi and Gamow- 
Teller transitions to (see §8.3), and via i/'e-electron scattering, labeled (solid line). The total number of Ue neutrinos detected via these channels 

is the thick solid line. The thick dashed line (labeled ‘All Others’) shows the contribution to the total neutrino signal from De, Vr, and Pr neutrinos via 

scattering on free electrons. For reference, we also show the Ve detection frequency {dN/dt, dotted line) during Ue breakout and the first 100 ms after bounce. 






